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Summary 
Dichloroethane-1.4 Dioxane at | atm 


Vapour-liquid equilibrium data has been determined for the two binary systems 


Monomethylamine-Dimethylamine at 10 atm. 
There is satisfactory agreement between the experimental activity coefficients and those calculated 
by the van LaAar equation 


Résumé—Détermination expérimentale du diagramme d’équilibre liquide-vapeur des deux binaires: 
dichloréthane 1.4 dioxane sous 1 atm 


monométhylamine-dimeéthylamine sous 10 atm 


Les résultats sont comparés avec ceux déduits de | équation théorique GIBBS DUHEM-VAN LAAR 


A. Systéme dichloréthane-1,4 dioxane 
I. Puretré pes CONSTITUANTS Les indices de réfraction ont été déterminés au 

Comme l’indique le Tableau | de constantes physiques, ™oyen du réfractométre Fery-V craBpoure donnant 

le dichlorethane et le dioxane utilisés ont été obtenus la précision de la 4°™° décimale. 

a un degré de pureté satisfaisant par rectification dans Les densités ont été déterminées au moyen d'un 

une colonne a grand pouvoir séparateur *. picnometre. 

Il. Tensions DE VAPEUR 


* La colonne dont nous nous servons couramment dans 
notre Laboratoire pour la purification de nos produits peut Le Tableau 2? reproduit les valeurs des tensions de 


étre considérée comme équivalente a 46 plateaux théoriques vapeur des constituants purs trouvées dans la littéra- 
s le cas ‘ re bine » veth yle yclohe rane 
lans le cas du mélange binaire + ptar ” yley rar ture. Pour interpolation nécessaire A nos calculs. 
avec un débit dans la colonne de 100 cc /h/em®*, débit maintenu 
constant sur toute la hauteur de la colonne dans des conditions nous nous sommes bornés a | emploi de la formule 
=> totes log P = f 7) traduisant par une droite sur un dia- 
Les analyses ont été effectuces apres une période de fone- 
de gramme A grande échelle. 


tionnement 8 heures. 


Table au ] 


Dichiorétane 


Diorane 


Experimental Littérature Experimental Littérature 


a 


réfraction 1.4445 1.4443 [1] 1,4226 1,4232 [1] 


Indice de 


Table au 2 


Constituant 1: Dichloréthane [4] 


Pression en mg Hg 10 20 40 60 100 200 400 760 


Température en °C ...... 44.5 24 13,6 2,4 + 10 + 18,1 + 29,4 + 45,7 + 64 83,5 


Constituant 2: 1,4 Dioxane [4] 


Pression en mmHg ...... l 5 10 20 40 60 100 200 400 760 


Températureen °C ...... 35,8 12,8 1,2 12 25,2 33,8 45,1 62,3 81,8 101,5 
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ANALYSE DES ECHANTILLONS 


L’analyse des échantillons a été faite par réfracto- 


métrie, la courbe n f(x) ayant été déterminée au 
préalable (Tableau 3 et Fig. la). 


Fig. la. Courbe nip 


Indice de 


moléculaire en 


20 
Tableau 3. refraction nj, en fonction du_ titre 


dic hlore thane 


In lice de 


refraction 


molec ula re 


T itre 
en dichloréthane 


0 1.4226 
1.4240 
1.4255 
1.4293 
14314 
1.4335 
1.4388 
1.4415 
l 1.4445 


0.100 
0.200 
0,400 
0.500 
0.600 
O.800 


IV. Mernope EXPERIMENTALE 

La détermination de la courbe d’équilibre liquide-va- 
peur a été faite en utilisant l'appareil de Stage [6] 
représenté sur la Fig. 2a. Le mélange liquide est 
chauffé a ébullition au moyen d'un chauffe ballon 
électrique recouvrant entiérement le bouilleur. Le 
distillat condensé est retenu par une cuvette annulaire 


entourant le col. C'est un appareil simple et robuste 


Etudes d’equilibres liquide-vapeur 


Chemical 
Engineering Science 


dont Vexactitude a parfois été contestée en lui attri- 
buant un léger pouvoir de fractionnement. Toutefois il 
a été constaté dans notre laboratoire au cours de 
étude du systéme eaudioxane que l'appareil de Stage 
a donné des résultats en accord avec la littérature [7]. 

La température est déterminée au moyen d'un 


thermocouple ferconstantan connecté & un potentio- 


Fig. 2a. stage. 


Appareil de 


métre portatif, capable de donner une précision de 
0.01 millivolt (0.2° C). 


faites au voisinage de 760 mm; mais le températures 


Toutes les experiences ont été 


du mélange en ébullition ont été corrigées et ramenées 
aux valeurs qu’elles auraient eues sous 760 mm. 
Une charge totale de 200 cc est utilisée pour chaque 
essai; le régime est établi au bout de deux heures 
d’ébullition environ. On préléve alors un échantillon 
de la phase liquide et un échantillon de la phase vapeur 


condensée (échantillons de 1 ce). 


V. 
Les résultats expérimentaux sont rassemblés dans le 
Tableau 4 et représentés sur les Fig. 3a et 4a. La 


Fig. 3a concerne le faisceau d’ébullition et de rosée: 
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Tableau 4. Données expcrimentales dé Uequilibre liquide-va pe ur du systeme dichloréthans -1,4 dioxane 


m perature °C ry V1 Veo log 


Vs log y 


100.1 0,05 0,130 1,445 l 0.15987 0 


98.7 0.110 0,220 1,285 0,995 0,10890 0,00218 
0,225 0.385 1,150 1,005 0.06070 0.00217 
4 0,325 0,495 1,105 1,008 0.04336 0.00346 
92.4 0.420 0.580 1,050 1.010 0.02119 0.00432 
0,555 O,715 1,040 0.01703 0,00432 
87.1 0.740 O.865 1015 1.025 0.00647 0.01072 
S4,8 0,900 0,945 l 1,065 0 0.02735 


Gu 


Fig. 3a. Courbes d’ébullition et de rosée systeme dichlor- 


éthane-1,4 dioxane. Fig. 4a. Courbe d équilibre liquide vapeur systeme 


dichloréthane-1,4 dioxane. 


la Fig. 4a traduit la relation y = f(x) indépendamment 


Les coefficients d’activité ont été caleulés d’aprés 
les relations suivantes: 
Py 
(La) 


Py 
P2 (1b) 


et les résultats obtenus sont représentés par les courbes 


de la Fig. 5a. 


VI. 


Les résultats expérimentaux ont été vérifiés au moyen y 


VERIFICATION ET CONTROLE 


af 
des équations de van Laar déduites de la relation ~<a 
générale de 


log 7 


log 72 


Bx, \2 ° Fig. 5a. Courbe log y:h(°) systéme dichloréthane- 
Ax, 1,4 dioxane. 
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ou A et B sont des constantes du systéme considéré 
dans un faible intervalle de température. 

Il convient de faire remarquer que ces équations, 
aussi bien celle de DuHEM que celles de vAN Laar ont 


été établies en partant de conditions isothermes et 


ne sappliquent strictement qu’a des résultats de ce 


type. Mais on sait qu'il est possible d’utiliser ces 
equations dans le cas de systémes dont les com posants 
ont, sous une pression determinée, des points d’ébulli- 
tion voisins, ce qui est le cas pour le systeme étudié 
ici, puisque lécart est inférieur A 20 

Les valeurs des constantes A et B ont été obtenues 


graphiquement en prolongeant les courbes représen- 
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tatives des coefficients d’activité en fonction du titre 


jusqu’aux points extremes ry et a, l. 
On trouve de cette facon 


B 0.050 


Les courbes résultant de lapplication des équa 
tions de VAN Laar A partir de valeurs précédentes de 
A et B sont représentées en pointillés sur la Fig. 5a 
L’accord entre ces derniéres et les courbes experi 
mentales (en trait plein) est aussi satisfaisant que le 


vrmet lemploi des équations (2a) et (2b). 


B. Systéme monométhylamine-diméthylamine 


I. Pureré pes ConstTirvants 


La monométhylamine est obtenue par dégazage d'une 
solution aqueuse 40% d'origine américaine. Toutes 
‘ les precautions ont été 


prises pour obtenir un 
produit pur et anhydre 
les vapeurs de mono 
méthylamine passant sur 
de la soude caustique en 
pastilles pour éliminer 
les entrainements d'eau 
éventuels. De plus le 
recepteur etait protege 
de Vhumidité atmos 
phérique par in plege 
contenant chlorure 


de ( ak ium 


If. Tenstons pe VaPeurR 
Nous rapportons dans le Tableau 5 les valeurs des 


tensions de vapeur des constituants purs. 


Ill. APPAREILLAGE 

L’appareil OTHMER nous ayant donné des mécomple- 
tes, nous avons préféré employer un appareil statique 
Il s'agit d'un autoclave (Fig. 1b) d'une capacité de 
1200 ce environ, muni d'un manométre et d'une gaine 
pyrométrique permettant le contrdle de la pression et 
de la température. Au moyen de deux vannes, on 
préléve simulanément d'une part un échantillon de la 
phase liquide par un tube plongeant, d’autre part un 
échantillon de la phase vapeur. 

On dabord le 


remplit ensuite jusqu’aé moitié (600 ce). 


fait vide dans l’'autoclave: on le 
L’autoclave 


est alors plongé dans un bain d’eau ou d’huile et 


Tensions de tpeur des constituants purs 


Pression en atm 


rempérature en 


Constituant 1: Monométh ylamine [8] 


106.3 133,7 144.6 


vapeur jous pre 


pour paise echant 
B Vai 
echantillom phase 


trique; D Tubs 


Vapour 
( Gaime pyvrom«e Pression en atm 


ur dans liquide; 


plonge 


Température en °C 
EE Manometr 


Constituant 2: Diméthylamine [9] 


132,2 162,6 164,5 


149.8 


La diméthylamine est obtenue dans les mémes 
conditions a partir d’une solution A 25%. 
La pureté des produits, déterminée par analyse 


chimique (acidimétrie) a atteint 99.9%. 


l'on chauffe jusqu’a atteindre la pression désirée. Le 
systéme est maintenu pendant trois heures environ a 
cette température et sous cette pression en agitant 


fréquemment (basculement) de fagon A obtenir un trés 
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bon équilibre des phases. On effectue alors une prise 
d’échantillon des deux phases, respectivement dans 
deux récipients contenant un poids d’eau connu et on 
dose lalcalinité de chaque échantillon. 

La vanne qui permet la prise d’échantillon de la 
phase vapeur est entourée d'une résistance chauffante, 
de facon a éviter toute condensation partielle. 

Le manométre, sensible A + 0,100 kg cm? est véri- 
fié et réétalonné avant et aprés chaque détermination 
Les températures sont enregistrées au moyen d'un 
couple fer-constantan relié A un potentiométre portatif 
donnant une précision de 0,01 millivolt (0,2 °C 

L’appareillage a été étalonné: 

a) en déterminant la courbe de tension de vapeur 
d'un corps connu, le méthanol 

b) en déterminant les données d’équilibre liquide- 
vapeur sous pression d'un systéme connu (systéme 
eau-méthanol sous 10 atm 
Comme en témoignent les Tableaux 6 et 7, les 


résultats expérimentaux sont en bon accord avec 


ceux données dans la littérature. 


Tableau 6. Tensions de vapeur du methanol 


rature 


T m pe rature 
Littérature) [1] 


Te m pe 


Pression absol ue 


en ka cm? observe 


Tableau A Equilibre liquide ra peur du syateme eau méthanol 


sous 10 atm. 


littérature [10] 


0,260 0.270 
0.440 
0,610 


0,150 


0.240 0.450 


O.615 


O45 0,740 0,740 


IV. ANALYSES 


Pour déterminer les concentrations dans l'une et 


l'autre phase on a procédé a des dosages alcalimétri- 


ques en se basant sur les équations suivantes: 


Si a, = poids du constituant 1 dans |’échantillon 
a, = poids du constituant 2 dans |’échantillon 
masse moléculaire du constituant 1 
masse moléculaire du constituant 2 

poids total de l’échantillon 
alealinité totale exprimée en grammes 


d’azote 


V. 


Les résultats expérimentaux sont rassemblés dans le 


Tableau 8 et représentes sur les Fig. 2b et 3b. La 


Courhe a a’é bullition et le rosee systeme mono 


Fig. 2b concerne le faisceau d’ébullition et de rosée; 


la Fig. 3b traduit la relation y f(x). 


Les coefficients d’activité ont été caleulés a l'aide 
des équations suivantes: 
Py, 


Py * 


Z, (3a) 


(3b) 
En effet, la vapeur ne suivant pas les lois des gaz 


parfaits, on doit corriger les équations (la) et (1 b) par 


un facteur Z qui tient compte de la compressibilité. 


Voll 
on 
P= a, + 4;, 
4 a 
M, M, * 
oe od’ l'on tire a, et dy. puis 
OL. 
952 
He 64,5 64,7 
3 
4 104,5 
5 113.4 113 
6 119 119 
» 
Ss 130.8 130.5 
10 138,7 138.5 
12 146 146 
15 155.9 155.5 
; 
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Tableau 8. Données expérimentales d'équilibre liquide-vapeur du systéme monométh ylamine-diméthylamine 


Température en °C log log 


0.070 0,130 1.260 ] 

1,125 1015 
1.055 1,038 
1.045 1.048 


1,025 1.105 


0.10037 0 


0.02325 


0,265 0,380 0.00647 


0.460 O.580) 0.01620 


0,620 0.02036 


0,730 0.01072 


0.04336 


big. 3b. Systeme monométhylamine-diméthylamine, 


D’aprés SCATCHARD et GILMANN [11], le facteur Z 
est calculable de la fagon suivante: 
P)(v — B) 
RT 


0 
InZ 


(4a) 


volume moléculaire du liquide 
second coefficient du viriel de léquation 
d'état, dépendant de la pression réduite et de 
la température réduite et défini par 

RT, 0,400 


B P 0,197 —0,012 


0,146 
he oa 
T; T, 3,27 
Ce facteur Z peut également étre lu sur un dia- 
gramme (12). 
Les valeurs de y, et y, ainsi calculées ont été rap- 
portées dans le Tableau 8 et l’on a tracé sur la Fig. 4b, 


les courbes traduisant la relation log y 


(2). 


VI. VERIFICATION ET CONTROLE 
De la méme facon que pour le systéme dichloréthane- 
dioxane, on a vérifié ces résultats expérimentaux au 


moyen des équations de vAN LaaR (2a) et (2b). 


Les valeurs de A et B, valeurs de log y, et de log y, 


aux points extrémes, déterminées graphiquement sont 


Fig. 4b. Systéme monométhylamine-diméthylamine. 
les suivantes 


0.140 


0.090. 


Les courbes résultant de application des équations 
de VAN Laar A partir de ces valeurs de A et B, sont 
représentées en pointillés sur la Fig. 4b. L’accord 


avec les courbes expérimentales est satis-faisant. 


CONCLUSIONS 


Le systéme dichloréthane-1,4 dioxane n'est pas un 
systéme idéal, non plus que le systéme monométhyl- 
amine-diméthylamine. 


Les légéres différences que l'on observe dans les 
courbes des coefficients d’activité entre les valeurs 
expérimentales et les valeurs caleulées d’aprés les 
équations de VAN LaaR sont dues en partie a l’ap- 
pareillage, et en partie A la non isothermicité des 
laccord est trés 


séries d’expériences. Néanmoins 


acceptable. 
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Summary. Experiments on the displacement of sugar liquors by water from columns of granular bone 


charcoal (sweetening off) have led to a theory which accounts for the channelling which sometimes oc: urs 


when one fluid follows another along a uniforml, packed column. The existence or absence of a tendency 


towards channelling is shown to depend upon the linear velocity of flow. A critical velox ity can be defined 


in terms of the viscosities and densities of the two fluids. The interface between the two fluids may be 


(1) inherently stable, (2) inherently unstable, (3) stable or unstable according to the relation between the 
actual and the critical velocities. In the case of “sweetening off’’—an important operation in sugar refin- 
ing—channelling occurs if the velocity of flow exceeds the critical ve locity. Experimental evidence is 


quoted in support of the theory. 


Résumé— Etudiant le déplacement par l'eau des jus sucrés dans les colonnes & noir animal, l’auteur propose 


une théorie qui rend « ompte des court-circuits qui se produisent parfois lorsqu’un fluide en déplace un autre 


dans une colonne a garnissage uniforme I] définit une v itesse < ritique fonction des viscosités et densités des 


deux fluides 


La surface de s« paration de deux fluides peut étre: 1, essentiellement stable 2. essentiellement in- 


stable 3. tantét stable, tantét instable, suivant le rapport entre la vitesse actuelle et la vitesse « ritique 


Dans le cas de la décoloration des jus sucrés, le court-cireuitage se produit si la vitesse d’écoulement dépasse 


la Vitesse ¢ ritique 


La théorie parait en accord avec les faits experimentaux cités par l’auteur. 


INTRODUCTION 


of market quality granulated sugar, the cistern is 


Refining of raw sugar, beet or cane, necessitates a se- Sweetened off—i.e. the sugar liquor remaining in the 
quence of processes which are devised to reduce the cistern is displaced downwards by water. 
concentration of non sugars. After affination the raw Any mixing of the liquor and the water during 


sugar is dissolved and the solution is purified by sweetening off is undesirable. The liquor to be dis- 


carbonatation and filtration. Before crystallisation placed usually contains about 67° of sucrose. Ideally 


on of granulated sugar by boiling in vacuum pans the the percentage concentration (~ Brix) of the effluent 
5 sugar liquors are further purified by downward per- should follow the broken curve of Fig. 1, the discontin- 
4 colation through granular bone charcoal (char) indeep uous reduction to zero Brix taking place when the 
- vertical cylindrical ‘“‘cisterns’’. In this process ad- water-liquor interface reaches the bottom of the 
: sorption of non-sugars, colouring matters, salts etc. cistern. 

takes place and as percolation proceeds the adsorptive In practice some departure from this ideal cannot 
: power of the char is progressively reduced and the _ be entirely prevented. The continuous curve is typical 


effluent from the cistern becomes less pure. When of practical cases. In the Tate and Lyle Refineries 


the purity of the effluent, as indicated by its colour it has been customary to accept liquor for crystallisa- 


density, has decreased to an established limit, beyond tion until the concentration falls to about 45° Bx. 


which the liquor would not be suitable for production Liquors below 45° Bx can be used for sweetening 
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off other cisterns [1]. Alternatively they can be 
concentrated by evaporation, in which case it is 
important to minimise the quantity of light liquors 
(liquors of concentration less than 45° Bx) in order 


to economise in time and cost of evaporation. As 


also the evaporated lights are of much lower purity 


Fig l Sweetening off Ideal and actual irves 


than the liquor from which they originated, through 
desorption from the char of previously adsorbed 
impurities, they cannot be utilised directly for pro- 


duction of granulated sugar. 


Propuction oF Licut Liquors DurING SWEETENING 
Orr 

Hitherto two main factors have been recognised as 
contributing to the formation of lights. These are (1) 
mixing of the water and the liquor by diffusion at the 
horizontal interface and (2) the time taken for the 
sugar liquor to drain and diffuse out of the pores of 
the granules. Turbulence has been mentioned as a 
third contributory factor [2] but we consider that this 
is improbable. The Reynolds number under Refinery 
conditions of operation is well below the limit at 
which turbulent flow can arise. As a matter of long 
established experience the quantity of lights can only 
be kept down to a reasonable level by sweetening off 
with a very small velocity of flow. This indicates 
that diffusion of liquor from the pores of the granules 
is the limiting factor. 

Recent laboratory and pilot scale experiments [3] 
have shown that the rate of adsorption of colour 
(impurities) during percolation can be greatly increased 
by using char of smaller granular size than has been 
usual hitherto. Therefore there is a considerable 
incentive to use fine grist char. Considerations of the 
pressure gradient required to generate the flow of 


liquor will limit the extent to which the char size can 
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be reduced. A limiting mean aperture (M.A.) of the 
order of 0-01 in. is imposed by this factor. However 
a more serious difficulty has been encountered during 
pilot scale experiments with fine grist char [4]. viz.: 
the excessive quantity of lights formed when sweeten- 


normal 


her 


Fig. 2. Sweetening off—Pilot plant curves for normal 


grist and fine grist char 


ing off proceeds at a rate of flow that is satisfactory for 
the normal grist of char. 

Two char cisterns were operated in parallel. One 
of these contained char of normal size and the other 


- 


Fig. 3. Sweetening on— Pilot plant curves for normal grist 


and fine grist char. 


a fine grist char. The size distributions of the chars 
were approximately Gaussian (normal) with respect 
to weight. The mean aperture of the normal char 
was 0-039 in. and that of the fine char was 0-011 in. 
The coefficients of variation (C.V.) were about 52% 
(normal) and 35% (fine). Figs.2 and 3 show. re- 
spectively, the curves obtained when the chars were 
sweetened off and sweetened on. Sweetening on is 
the inverse of sweetening off. Char may sometimes 


be settled into water. Sweetening on is the displace- 


a 7 
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ment of this water by sugar liquor, in this case by 
downward flow. 

We see that at the flow velocities used the normal 
char gives a satisfactory sweetening off curve and an 
unsatisfactory curve for sweetening on, while the 
opposite is the case for the fine char. A theoretical 
explanation of these phenomena is advanced in the 
following. 

THEORY OF CHANNELLIN, 
When a mobile fluid (¢.g. water) displaces a viscous 
fluid (¢.g. sugar liquor) from the voids of a uniform 
packing the difference between the viscosities tends 


to cause instability. 


Fig. 4 


Consider the interface A BC D E (Fig. 4) between 
the two fluids. BC D is a bulge which is supposed 
to have developed in the interface. The line GC F 
which is perpendicular to A E is a path of minimum 
resistance to flow so that the velocity through C ine- 
reases and the depth of the bulge also increases. There- 
fore the difference between the viscosities tends to 
make the interface unstable during the process of 
sweetening off. Conversely this factor has a stabilising 
influence on the interface during sweetening on. 

If, during the operation of sweetening off or 
sweetening on, the less dense fluid is above the denser 
one the difference between the densities of the fluids 
tends to keep the interface horizontal. On the other 
hand it tends to cause instability of the interface if 
the denser fluid is uppermost. 

Sugar solutions are denser and more viscous than 
water, and from a consideration of these two factors 
(density and viscosity) it follows that upward sweeten- 
ing off is an inherently unstable process while upward 
sweetening on is inherently stable. The other two 
processes— sweetening off and on by downward flow 
may be stable or unstable according to which of the 
two factors predominates. 

The theory may be formulated quantitatively as 
follows. 

Consider a fluid of viscosity mw. and density 0. 
displacing a fluid of of viscosity a, and density 0, 


by downward flow (Fig. 5). 
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Channelling packed columns 


The velocity of flow = r. 

Suppose that a shallow depression C D, of depth 
Ox, develops in the interface A B. 

Let py — pressure at the level A B, 
P pressure just below C D, 


P. = pressure just above C D. 


The depression will disappear if p, > p,, and will 
increase in depth if p, > p,. 
Now 
Py Po ~ 90,04 
and 


Ps Po k Mgt Ox 


where & resistance constant of the char. and g 
acceleration due to gravity 


Fig. 5 
The condition for stable flow is therefore : 


Py Py = 9(0, — 0,)0x— ke (My — Mg) Ox 0, 
That is 


Oo) k Me) 


or 
— G2) 
— fy) 
Apart from & all of the quantities in this expression 
are known. k must be determined experimentally. 
q(o, iy) 


will be called the critical velocity, 


If o, 0, and 4, > wy (downward sweetening off) 
the flow is unstable if » exceeds the critical value. 
In the case of downward “sweetening on” (0. > 0, 
and s# > sy) the flow is stable if (and only if) the 
velocity exceeds the critical value. 

Note that conditions that vield stable sweetening 
off cause unstable sweetening on and vice versa. 
Note also that laboratory experiments designed to 
test sweetening off and on must be done at the correct 
linear velocity. 

The differences between the curves for fine and 
normal char in Figs. 2 and 3 can now be accounted for 
by the difference between the two values of k. The 
fine char has a greater value of k and therefore a 
smaller critical velocity. 

Owing to lack of exact data the critical velocities 
could only be roughly estimated but Table 1 shows 
that within the limits of uncertainty the curves conform 


to the theory. 


Voll 
A 
— | 
A 8 
Po 
dx 
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*low | 
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Two direct tests of the theory have been made. 
In the first test a transparent perspex cell of height 

13 in.., 
containing sugar solution was packed to a depth of 


breadth 6in. and thickness lin. and 
9 in. with glass spheres of 0-4 mm diameter (Plates | 
and 2). The value of k for the packing was determined 
and the critical velocity calculated. The cell was then 
sweetened off by adding coloured water at the top 


™~ 


Fig. 6. Sweetening off—Curves obtained with the 1 ft 


laboratory column 


from a constant level device and running off the sugar 
the The cell 
illuminated and photographs were taken when the 
volumes of liquor displaced were |, }, j, 1 and 1 } times 


Plate | The dark 


areas of the cell are occupied by the coloured water 


solution from bottom was suitably 


the void volumes of the packing 


and the light areas by the sugar solution. Sets 
Table ] 
Sweeten- Liquor velocity (it hr) Result 
(irvat 
na 
Ietual Critical Expected Ohserved 
off (26 good yood 
normal 
on Ol bad bad 
off 0-26 bet ween ? bad 
fine 0-12 and 
on good good 
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numbered 1, 2, 3 and 4 of Plate | show the results 
obtained when the velocities of flow were respectively 
1, 1, 2 and 4 times the critical velocity. Satisfactory 


separation of the liquid is maintained in sets | 


(v and 2 but channelling occurs 
in sets 3 (¢ 2r.) and 4 (rv 4v.). Similar sets 
of photographs, shown in Plate 2, were made while 
sweetening on. In this case the results of sets 5 


with the I ft 


Fig. 7 Sweetening on—Curves obtained 
laboratory column 
Sv.) and 6 (¢ r,) are unsatisfactory, while 


sets 7 = and 8 (v 4¢,) show good separation. 


a column | ft high and 24 in. 


In the second test 
in diameter was packed with char of M.A 
and C.V 25%. The 


on at velocities 


and 
For 


convenience the experiment was made with 57° Bx 


char was sweetened off 


varving from tr, to 


sugar solution at 30°C. Under these conditions 


vr, -2 ft hr open tube velocity. The two yroups 
of curves, Figs. 6 and 7 show the results, which also 
conform to the theory. Satisfactory sweetening off 
is achieved only if ¢ < », and satisfactory sweetening 
on if 

The effects described are in dependentof small in- 


homogeneities in the packing. 


REFINERY CONSIDERATIONS 
We may now consider the quantitative nature of 
this problem as it affects Refinery operations. A char 


cistern in which liquor is decolourised by char will 


a 
| 
: 
q 
| 
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deliver a quantity of liquor which depends on the Pilot Plant Char Department at the Tate & Lyk 


time taken to pass one cistern volume of liquor. Ii Plaistow Refinery and to the Directors of Tate & Lyl 


this time is 3 hrs in the case of carbonatated liquor for permission to publish this paper incorporating 


on char of M.A. = 0-014 in. the volume of acceptable work carried out in the Research Laboratory. 


liquor produced is of the order of 20 cistern volumes 


NOTATION 


so that the cistern is decolourising liquor for about 


Bx % concentration of a sugar solution. 
1 6Ohrs. In the Refinery 67° Bx liquor is processed 
. M.A, Mean aperture Sieve aperture through which 
70°C. Such li has a critics tube veloci 
; at 70° C. Such liquor has a critical open tube velocity 50% of the char can be passed. 


through a packing of this char of about 0-4 ft/hr. C.V. = Coefficient of variation 
Assuming 60° voids in the char this gives a critical 100 « (Standard deviation) , 


linear velocity of 0-67 ft/hr. A separate experiment has M. A.) 
i Viscos of the fluid being displaced. 
Jensity of the flux elng displaced, 
does not fall to the value 1-8° Bx (at which it is : eye oe 
My Viscosity of the displacing fluid. 


turned to drain) until about 1-5 cistern void volumes 


Me Density of the displacing fluid. 


have been displaced. Thus the time taken to sweeten 


) Hydrostatic pressure at the plane part of the 
off at the critical velocity is 1-5/0-67 2-25 hrs/ft of interface. 
column. Py H ydrostati« pressure just below the bottom of the 
Conventional char cisterns are from 16 to 20 ft bulge. 
OL. deep so that sweetening off would take from 36 Po Hydrostatic pressure just above the bottom of 
I to 45 hrs. Comparing this with the 60 hrs of decolouri the bulge. 
sation we see that the non-effective time is a substantial g ~ Acceleration due to gravity. 
Resists ce sti ) MACK 4 *ressure 
fraction of the whole. On a shorter cistern the fraction Resistance constant of the packing Pressure 
, gradient required to generate flow at unit velocity 
would, of course, be correspondingly smaller , =? 
of a fluid of unit viscosity. 
rhe practical sweetening off of fine char is there- ; 
da Depth of an incipient bulge in the interface. 


fore still an incompletely solved problem. It is also 


Open tube velox ity of flow. 


one of considerable importance. The procedure 


Critical velocity of flow. 


customary with normal char has been found inad- 
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equate. Slow sweetening off of short columns may 


be an acceptable compromise but a really satisfactory 


solution is still awaited. 
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Equations based on theoretical considerations have been developed to give the volume trans 


Summary 
ported per minute, the loading in the drum, the time of passage and the power requirements for the trans- 


port of granulated solids through horizontal rotary cylinders without flights or any other lifting devices 


The simple equations are suitable for practical purposes. Experiments have shown that the results of the 


calculations are in reasonable accordance with the values obtained from the experiments. 


: Résumé Des considérations théoriques ont conduit A des equations permettant de caleuler le volume 

ie transporté par minute, la charge du cylindre, la durée de passage et la puissance consommeée pour le déplace 

a ment de solides en grains au travers de cylindres rotatifs horizontaux dépourvus de chicanes ou de tout autre 

ti dispositif élévatoire. Des équations simples suffisent pour les applications pratiques. Les résultats expéri 

in mentaux sont en bon accord avec les valeurs calculées 

1. INTRODUCTION calcining processes, in chemical reactions with gases, 
For handling a continuous flow of granulated material etc. rotary cylinders are often used, through which 


in drying, humidifying, mixing, heating, cooling or the granules are transported. 
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In many cases the interior of the rotary cylinders 
contains flights or other lifting devices causing an 


intensive contact between the granulated material 

and the gases that flow through the drum in counter- 

current or in parallel current with the granules 
The most simple construction of rotary drums 


however, is that without any flights or lifter blades 


inside the cylinder. This type of cylinder is often used 
in industry, e.g. for the manufacture of cement, for 
burning light granules of clay and so on 
The design of rotary cylinders involves the calen- 
lation of 
SECTION 6-8 


Surtace 


of material flowing through the 


1. The 


drum: 


amount 


2. The loading in the drum as a percentage of 
the total volume: 

3. Time of passage of the material through the 
drum 

These quantities have to be determined as functions 
of 

1. Dimensions and other constructional data of the 
drum 


transported 


2 Properties of the material to be 
3. Operating conditions of the plant 


In horizontal drums, without flights and without 
flow of gases transport of granulated materials takes 
place because of the difference of the levels of the 
material at the entrance and discharge end of the 
drum [1]. Caleulation of the volume transported pet 


minute ®, gave the equation 


h, DAn.cot 
(A) 
“\D L 


where a and } are constants. For the meaning of the 
other symbols, see list of notations. 

This result contradicts the empirical equation of 
SULLIVAN, et al. [2] according to which the transport 


would be nil in horizontal rotary kilns. 
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Last year SAEMAN published a theory of the trans 
port of material through rotary cylinders [3]. For 
lightly loaded kilns an equation is derived resembling 
that of Suttivan. Due to the simplifications intro 


duced into the calculations, theoretical 
give results than can be used 
For the 


of material at heavy kiln loads and in horizontal 


considerations do not 
for horizontal rotary cylinders. transport 
evlinders no solution is given 
Eq. (A) gives figures for the transport of material 
through horizontal rotary cylinders, that are in good 
agreement with the experiments, even for heavy 
SECTION A-A loads. It is interesting to note that 
this equation has been developed 
on the same theoretical assumptions 
as used two years later by SAEMAN 
The different result is mainly due to 
the fact that the simplifications 
necessary for the mathematical 
development of eq. (A), have been 
chosen more fortunately than those 
The 


calculation of 


used by SAEMAN. results of 


the theory for the 
material transport, load, time of passage and power 
1949 


1}--can immediately be used for practical purposes 


requirements,—already published in part in 


and will, therefore, be given in some detail in the 


present paper. 
2. Egvation vor Transport Rate oF Marertal 
The surface of the granulated material in a horizontal 
rotary drum without flights is shown schematically 
in Fig. 1. The feed rate should be sufficient to keep 
the bed depth at the entrance constant (measured 
along B-B) 


slide along the cylinder wall 


It is assumed that the particles do not 


The particles cascade quickly down ina thin layer 
on the surface of the bed along chord k and are then 


carried back again in the bulk of the bed to the surface 


The position of the surface with respect to the hori 


zontal plane is determined by the angle of repose / 
of the material. The angle of repose depends on the 
physical properties of the material, as friction coef 
ficient of the particles, distribution of size of particles 
and bulk density. The axial component of the path 
of the cascading particle can be determined by geo- 
metrical considerations as a function of the path k 
the angle of repose § and the angle e: 
_ cos p.sin tane 


tan* + sin? 


ié 
3 
a 
F 
2 
- | 1952 
a 
Fig. 1 of moatenal i horizont tary cy 
fe 
4 
4 
294 
= 
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The angle ¢ being small, tan* ¢ may neglected 


against sin® so 


cot p . tan é 
where 
dh 


al 
tané da 


A particle on a circle with radius ~ travels along the 


are gr in the time When the particle arrives 


at the surface of the bed, it cascades along a path of 
a length equal to chord & in a time that may be neg- 
lected with respect to the time needed for travelling 


along the are of the circle. The average axial transport 


velocity v, of this particle during the complete cycle 


will therefore be 
dh 
A da 


cot p 


a CTOss- 


This equation is valid for all particles in 


section of the area qr -d and the flow of material 


through this elementary section is 


dh 
da 


n.cotp. (3) 


AD. 


The total flow of material through a section of the 


drum is 


dh 
2an.cotp. Jk. dr, (4) 
dz. 
where (R —h)*. 


Integration of (4) gives the equation 


4 h h \2 R 
7 R*n cot R | (5) 

Writing and > where the sub 


scripts indicate values at the entrance and the end 


of the drum respectively, the eq. (5) can be inte 
grated 
L 
| 2&— —F(§,). (6) 
ia R*.n. cot ak ‘ 
0 


In cases where h,; — 0 the function F(g,) becomes: 


| 2] 4 (7) 
(25, so 2>0 So) | 
Fig. 2 shows F(&,) as a function of &. It can be seen 
in Fig. 2 that the rather complicated expression (7) 


can be replaced by the sufficiently close approxi- 


mation : 
£2-08 


F (&,) = 0-62 . 
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in the interval 0-3 0-9 that is interesting for 
practical purposes. 
Substitution of this approximation in eq. (6) 
gives for the volume transported per minute 
2 OR 4A n. cots 
— 0-68 ("° (8) 


for values in the interval 0-15 


0.4 


208 
—— CURVE FOR 0,62}, 
@ VALVES CALCULATED FROM EQ(7) 
0,02 


+ 


001 
01 02 os 
Fig. 2. Curve representing function F'(&, 
3. Laporatory EXPERIMENTS 


a) Description of apparatus 
For the experiments, glass and steel cylinders have 
both of 


The constriction at 


been used, with a diameter 110mm and 
a length of 650mm (see Fig. 3). 
the entrance was variable, and so experiments could be 
carried out with different values of h,. The discharge 
end of the rotating cylinder was not constricted and 
no experiments have been made with such constric- 
tions. A constant supply, necessary for the experi- 
ments, has been obtained by a simple but reliable 


control device. 


of here 


fication proposed by Prof. H. Kramers of the original de- 


* The development eq. (3) given is a simpli 


velopment by the authors. 
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b) Qualitative ea pe riments 

Two materials have been used for the experiments 

riz. dry sharp sand (diameter of particles 0-2 to | mm) 

and rve (particles 7mm long and 2 mm diameter) 
For facilitating the study of the cascading 

of the granules the first experiments have 

been carried out in a glass cylinder. It turned 

out that in the glass cylinder the bed of sand 

and even more so the rye, started sliding as 

one solid mass as soon as the depth h ot the bed 

(h 


took 


became less than a certain value 12 mm 


and no cascading of particles place 


used for the experiments 


because the angle of repose of the material was 


not attained. 
In the steel evlinder with a rough surface the 
particles cascaded, without sliding of the material as 


16 


— 


——SAND, 2 = 149kg/! 
{20,71 


nN 


- (LITRES 


be 


a nirpm) 


Fig. 4. Volume flow of transported materia! 


a whole, because the friction between the material 
and the wall was greater than the mutual friction 


between the granules. 

The movement of the material in a cross section 
could be observed by using a transparent ring of 
At 


celluloid at the entrance of the cylinder. low 
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speeds of the cylinder a particular granule described 


continuously about the same path. The layer of cas 
eading granules gave the impression of a laminar 
flowing liquid. At higher speeds, x > 15 r.p.m., the 
straight cascading path changed 
into a wave line with a minimum 


in the middle. 


¢) Quantitative experimne nts 
When 
evlinder and a certain depth /, 
ot hed 


heen set up, the feed of mat rial 


a certain speed of the 


the at the entrance had 


was controlled in such a way that 
ismall amount of material flowed 
over the ring out of the cylinder 
state was ob 


The steady soon 


iined and then the transported 
collected for 
The 


accuracy of the measurement was 


material would be 


some time and weighed. 


2 
Fig. 4 


measured 


ibout 


of the 
quantities trans 


The load of 


the rotating cylinder was determined by weighing its 


shows curves 
of 
ported material. 


contents at the end of each experiment 


d) Bulk den ity 


In the beginning of this study the bulk density has 
heen determined by means of a calibrated measuring 
glass. The values obtained in this way showed large 
variations, depending on the manner of filling the 
measuring glass. Careful filling or prolonged stamping 
of the contents resulted for sand in values varying 
1-45 to 1-62 


measuring the bulk density has been adopted. 


Therefore another method of 
A eyl- 
glass discs at 
of 


After rotating the cylinder for some time 


from 


inder of 110mm diameter, closed by 


both ends, was charged with a certain amount 
material. 
the depth of the bed of material was measured accu- 
rately and the volume of the bed calculated. From 
the weight of the contents of the cylinder the bulk 
density under actual conditions in a rotary drum 
could then be computed, 

The bulk density measured in this way was for 
sand 1-49 kg/l and for rye 0-71 kg/L. 

The angle of repose § has been measured in rest. 
After slowly revolving the cylinder by hand until 


cascading of the granulated material started, the 
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and a _ horizontal 


bed 


plane could be determined with an accuracy of | 


angle between the surface 


The value of the angle of repose was for sand 33 


and for rye 40°. 


4. COMPARISON OF THE RESULTS OF CALCULATIONS 
AND EXPERIMENTS 
The results of the experiments are given in Fig. 4, 
which shows for each value of /,/ D a linear relation- 
ship between the volume transported per minute and 
the speed of the rotating cylinder. From these results 
an experimental coefficient /,... can be determined for 
each value of h,/D (see Fig. 5) 
L hy, 


n.cotp 


Comparing the results of the experiments, shown 
in Fig. 4, 


(eq. (8)], it appears that the experimental values are 


with those of the theoretical calculation 


approximately 13° higher than the theoretical ones. 
As the measured values in Fig. 5 are lying rather 


close on a straight line, the following empirical 


equation can be given for the volume ®, transported 
per minute 
2-05 


DP, ‘ D L 


h, 
In 


eq. (S), except for the exact values 


which equation is for values 0-15 


accordance with 
of the coefficients. 


DruM 


bed 


CALCULATION OF THE LOADING IN THI 


cross-section of the 


be area of 


If X the 


(Fig. 1), then the load 7 in m® of the rotary cylinder is 


any 


H=—(X. (10) 
The calculation of H is simplified by using an ap- 
proximation instead of the exact formula for the area 


of the segment X. The equation 
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The load ratio of the rotary drum is defined as 
loading in cylinder 

where X, is the segmental area of the cross-section 

through the bed at « = 0. From (11) and (13) the 

value of 7 can be calculated 


0-67 LP th, 1.42 L 


O-5Y 
114227 (h,/ L 
0,2 

0.10 

0,08 

0,06 
004 
w 

f= 

‘ O*n cotg? 

© EXPERIMENTAL VALVES OF (SAND) 
x EXPERIMENTAL VALWES oF (RYE) | 
0.01 
0.10 02 04 06 O8 10 
Fi ». Comparison of calculated and measured 


coetiicients 


This shows that the load ratio has a constant value 
independent of size and operating conditions of the 
rotary cylinder and independent of the properties of 
the transported material. The loading H in the cevlinder 
has been determined in some of the experiments. 
Table | 


the values of 7 as calculated from the measured data. 


shows the results of the measurements and 


i 
X ~ 1-14 D?| j (ll Table I. Expe rimental values of the load ratio 
} D 0-20 0-29 0-38 
gives for values A/D < 0-45 an approxima . 
tion with an accuracy of 3°. 
ry 7] 6 Ss -p. 
The flow of material @, being the same ae 0 a 
Hex, 0-52 0-55 1-00 1-07 1-47 1-55 litres 
through every cross-section it can be found 
X,.L 0-90 0-90 1-55 1-55 2-30 2-30 litres 
. Nexp 0-58 0-61 0-645 0-695 0-64 0-675 
h \2-05 L—a 
(12) 
ho L 


By substituting eqs. (11) and (12) in eq. (10), the 
latter can be integrated: 
142 


° hy . 
H = 0-67 D?| (13) 
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6. CALCULATION OF TIME OF PASSAGE 
The average time of passage T of the particles is 


H 


D tan 
T 0-9] 


(min). (14) 


DP, dD hy 


e 
7 
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7. CALCULATION OF POWER REQUIREMENTS 
A particle cascading along a chord of length & (see 
cross section, Fig. 6) is lifted again over a height of 
ksin 8. The amount of work done on a particle 
entering the bed at P and leaving it at ( is 


y.r.dg.dr.dx.k.sinp (kgm). 


Per revolution, the volume entering at P is 
dr.dx f 


The power required for displacing the material in 
a bed segment of thickness dx, at » r.p.m. is found 


by integration 


dp d 


go sinp. 


da 
R-h 


material 


Cross-section thre ugh bed ot 


The solution of the integral in this equation is 


With an accuracy of 5% this integral can be replaced, 


for values 0-2 AVR < O8, by the equation 
2 h A \2 4 h 
Using this approximation the power requirements 
are 
an 2 
p 20) y.sinp R? ‘24 hdx (15) 
after substituting eq. (12) the integration of (15) 
gives 
(metrich.p). (16) 


S. CONCLUSIONS 


a) Through a horizontal rotary drum with a depth 
h, of the bed of material at the entrance and with a 
free discharge at the other end of the drum. material 
is transported as a result of the difference of the 
levels of the material at both ends. 

b) From theoretical considerations an equation for 
the nett volume transported per minute has been 


developed for values of A/D in the interval 0-15 
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< 0-45: 


. cot 
0-68 | 


L 


Experiments gave about 13% higher values: 


(m*/min). 


2-05 “wot A 


4 
ex D L (m*)mun). 


c) The load ratio 
H 

is 60 to 65% and is practically independent of the size 
of the drum, the operating conditions and the prop- 
erties of the transported material. 

d) The time of passage of the material through the 
drum is given by eq. (14) 


tans 
T= 09] 


D/ 


(min). 


e) The power needed for lifting the particles in 
the bed in the drum up to the cascading laver is 
given by eq. (16) 


p=O19n.y. L hg. sing (metric h.p.). 


NOTATION 
axial component of the path of a cascading particle 
angle between drum axis and surface of material 
(Fig. 1) 

D angle of repose of material 

bulk density of material (kg /m*) 

h depth of bed 
ho depth of bed at entrance of drum, 


L length of drum (m) 


r axial distance of section to entrance of drum 

X cross-sectional area of bed 

X, = cross-sectional area of bed at entrance (x — 0) 

i length of path of cascading particle along chord of 

radius 7 

q angle on chord k 

R D2 radius of drum (m) 

} speed of drum (revolutions per minute), 

v, average axial transport velox ity of a single particle 
Pp, volume transported per minute (m*® min) 


H = load of drum (m!*) 

y load ratio of drum = H/X,. L 

T= time of passage of a particle through the drum 
(min) 


p = power (metric h.p.) 
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2 The passage of granular solids through inclined rotary kilns 


Receive l 


Summary 
as first indicated by VAHL, 


end constrictions 


Résumé 


The passage of granular solids through inelined rotary kilns 
H. Kramers and P. Croockewit 


Laboratorium voor Physische Technologie, Technische Hogeschool, Delft 


Auqust 


En-utilisant la théorie simplifiée proposée par Vani pour le déplacement des solides dans un 


1952) 


Starting from a simple mechanism of the movement of granular solids in a rotating cylinder, 
formulae are derived for the hold-up of an inclined rotary kiln, with and without 


A good agreement is found between the theoretical values and the « xperimental results. 


cylindre tournant, auteur établit des formules rendant compte de la rétention dans un four rotatif incliné, 


muni ou non de chicanes. 


In the design of rotary kilns one of the main problems 


INTRODUCTION 


is the calculation of the kiln dimensions and operating 
conditions for a given feed rate, and certain require- 
ments as to the time of residence of the material. 

The first attempts to find a correlation between 
these variables were made by SuLiivan, Marer and 
taLsToN [1] who experimented with small inclined 
kilns and relatively small feed rates. They derived 
an empirical formula for the time of passage through 
kilns without constrictions, whereas several correc- 
tions were indicated for use in kilns with constrictions. 
This well-known work has been followed by minor 
adjustments by Bayarp [2]. Recently, SAEMAN [3] 
published a paper in which a sound theoretical basis 
was put forward for the calculation of kiln hold-up, al- 
though he did not succeed in giving a general solution 
for practical use. 

The theoretical results obtained in the present 
work are based on essentially the same mechanism 
of transportation as put forward by Sarman. They 
must, however, be considered as an extension of the 
ideas of Prof. VAuL, the result of which was pub- 
lished in 1949 [4]. A full account of VAHL’s work for the 
case of horizontal rotating cylinders appears in this 
journal [5]. In our derivation of the hold-up relation- 
ships for inclined kilns we will refer as much as possible 


to the latter paper. 


2. Derivation or FoRMULAE 

Transport mechanism 

We start from the schematized description of the path 
of a particle which is transported through a rotating 
cylinder, as proposed by VAut [5] and by Saeman [3] 
(see Fig. 1): The bulk of material in the kiln has no 
relative motion with respect to it. So a particle at 
a distance r from the kiln axis describes a circular 
path in a plane perpendicular to the axis with the 


same angular velocity as the kiln wall. After it has 


Les valeurs expérimentales sont en bon accord avec celles déduites de la théorie. 


reached the surface of the material it cascades down 
along a straight path having a “dynamic angle ot 
repose "’ with respect to a horizontal plane. The latter 
movement has two components, one perpendicular 
Only 
makes a contribution to the transport through the 


kiln. 


The ratio between s and k depends on the angle 


to the kiln axis (4) and one parallel to it (s). 


of repose (/) of the material, the slope of the kiln 


- 
A 

Za" 
i- 
mA 
g 
i 
Fig. | Schematic representation of rotary kiln with 


svimbols used 


(tan x) and the gradient of the surface with respect 
dh 
dx 
found from geometrical considerations and has been 


to the kiln axis This relationship can be 
shown by SAEMAN to be: 


cos dh 


k sin sinf 


tan % 


(1) 


It is valid for the small values of « encountered in 
dh 
dx 


With the same reasoning as used by VAHL ([5], 


practice and for small values of 

eq. (2), (3) and (4)) we arrive at the following formula 

for the volumetric flow of material through any cross- 

section of the kiln: 

tans dh h 

- ‘ 


3 sin p 


Pp, 


which is the basis for further calculations. It corre- 


sponds with eq. (18) in SAEMAN’s paper. 


Height of material along the kiln 


Eq. (2) is a differential equation which relates the 


height of material A with distance x from the feed end 
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for given conditions of operation As its solution is 
rather complicated we have simplified eq 
replacing the function containing h R by a linear one 

h A? 


» 
? 1-24 


- 


Fig. 2 


illustrates the degree of approximation ob- 


Fig. 2 
tained. The introduction of this approximation before 


Fig. 3. Dimensionless representation of the relat 
material and distance trom d 


tween height of 


end, for various values of 4; 


integrating eq. (2) is at variance with VAnt’s deri- 
vation for horizontal kilns (« 0). 


Further. if we use the following dimensionless 


combinations: 
R« op 


and N, (4) 


sin B 
No 


n tan x 


e of granular solids through inclined rotary kilns 
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it can readily be seen that eq. (2) contains two Vari- 
ables only, vez. 
h 
R Ne an L Np Av 
The solution of (2) is fixed, if the terminal condition 


for the height of material is given 
h hy, for r L. (5) 


We assume h, to be equal to the height of an eventual 
diaphragm constriction. If no constriction is present 
at the discharge side of the kiln, hy = 0. 

After substitution of (3) and (4) and for the con- 


dition (5), eq. (2) gives the s ution 


0.1903 
R N iD L 
j 


hy 
RN, RNo 


0-193 In (6) 


0.193 
RNg 


This rather complex relationship is represented by 
Fig. 3. It 


which the material height 4 changes with increasing 


shows, in dimensionless form, the way in 


distance from the discharge end, for various values 
of 

A special case presents itself when A, RNg 
0-193, giving a constant height over the whole length 
of the kiln. 
(low feed rate and long kiln with great slope) A, ap- 
proaches the value of 0-193 RN4,. 


the conditions at the discharge end represent a rela- 


For sufficiently small values of N, \, 


For such cases 


tively unimportant end effect. 


Calculation of fractional holdup X, 

Once the relationship between / and « has been 
established, the calculation of hold-up involves no 
difficulties 
filled with material to a uniform height /, the frac- 


X, of the 


material is 


essential If a cylinder of radius R is 


tion cross sectional area occupied by 


whereas 


This again is a complicated relationship between X, 


and h R. So we approximated it by the formula: 


(5) 


Fig. 2 shows the significance of this approximation. 
Actually a closer fit would have been obtained by 
a slight modification of the constants in (8) (compare 


[5], eq. (11)). However, we preferred the approxi- 


“it. 
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because of its advantage for 


here 


mation given 
numerical calculations. 
Since / varies through the kiln, the total relative 
hold-up X, has to be found by integration 
ée 9 
0 
was carried out with 


This operation graphically 


eq. (6) as a starting point. Because of the approxi- 
mation (8) the results are obtained in the form 

hy 


As we preferred to present a relationship between 


{(No N,) for different values of 


the two main variables X, and @,. we rearranged it 
so as to give a relation of the form 
X,N, (N, N,) 
hy 
RN, 
Table 1 gives the result of these calculations. 


(10) 


for different values of 


Table J. Values ol 


The passage of granular solids through inclined rotary kilns 


OOS 


0-10 


O21 0-021 0-022 0.022 

] 050 0-050 0-053 0-056 0-061 0-082 

O110 O117 0-127 06-140 0-168 0-254 

5 0-27 O29 O33 O40 O50 O61 1-07 

O51 O58 O71 O95 127 2-90 
20 O92 106 1158 2-28 413 7-65 


The relationship (10) is represented in Fig. 5 by 


the solid line for hy; = 0, i.¢. for a kiln without a con- 
striction at the end. For low and for high values of 
N, the curve approaches to different asymptotes. 
The former one corresponds to the case where the 
discharge conditions of the kiln are negligible, as 
For that 
be caleulated directly from (8), X, 


X, can 


X, and RN, 


mentioned before. case therefore, 
0-193, giving: 


X ODS (Ny N, jis for (ll ) 


No N, < O1, 
The asymptote for high values of N, N, refers 
to the case where the axial slope of the surface 
(dh/dx) is mainly responsible for the transport, e.g. 
when the kiln is nearly horizontal. Neglecting the 
term containing « in eq. (2) and using the proposed 
approximations, one finds: 
A, N; 


0-158 (Vo N, 23 for 


No oO 


and =0| 


(12) 
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or X, 


158 (Ny 3 for and hy 


The latter equation no longer contains x as a variable. 


Effect of end constrictions on hold-up 

From Table | it can be seen that for relatively small 
values of NV, N, the effect of an end constriction is 
only small, because for all values of h,, the eq. (11) is 
then approached 
effect of a 
a correction on the total hold-up calculated for 
hy, = 0. 


asymptotically. In general, the 


constriction can be regarded as being 


A correction factor f,, was defined by 


(X,),, (X 13) 


Fig. 4 gives a plot of /, as a function of N, \, for 


different values of h, RN,. Thus the total hold-up 


ojo 
NaN 
Fig. 4. Chart for evaluating the effect of an end constri 
tion on the hold up as found from Fig. 5 (hy 0 


of a kiln with end constriction. or with a given height 


he ab the dicchares end. can be calculated by first 


deriving (X,),, 9 from Fig. 5 and afterwards correct- 
ing it by using Fig. 4. Both graphs are based on the 


numerical results contained in Table 1. 


3. EXPERIMENTAL Parr 
After the above calculations had been set up tenta- 
tively, we compared them with Suttivan’s data for 
hy, = 0 (Tables 1, 2, 5, 6, 7 and 8 in reference [1}). 


The agreement appeared to be very good, as can be 


seen from Fig. 5 (crosses). It was then decided to 
carry out some more hold-up measurements, in 


particular in the range 2 No 20, i.e. at 
relatively vreater 


eft al. 


loads than those investigated by 
SULLIVAN 

The experiments were performed in a rotating 
0-197 m and 


the 


cylinder with an inner diameter 2R 
length L = 178m. It 


of a semi-technical 


was mounted on base 


cement kiln having the means 


for adjusting the kiln slope « between 0-5 and 5 
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and the rotational speed n between Table 3. Experimental results for hy 0 


0-25 r.p.s.*. The necessary 
0-04 and 0-25 r.p.s.*. The necessary Slop n Feed rate Hold-up Ni? (X 
precautions were taken to ensure tan (r.p.s.) (107% kg/sec) (ka) No Ne (Xp)eate. 


a constant feed rate of material. In 


each experiment, after steady con- Sand 
ditions had been reached, the feed 0-059 5-15 8-10 0 101 0-795 18-3 O-111 
was stopped and the hold-up of the | 0-080 2.68 5.00 0-062 (0-488 6-25 0-054 
= cylinder was determined by weighing odes 0-195 13-2 775 0097 0-763 14-2 0-088 
. its contents. We used river sand | (0-232 7.2 3.25 0048 0-378 6-55 0-043 
. and ground coke. The properties of 
0-040 6-38 13-3 0-166 1-22 29-7 0-169 
= 0-040 5-00 11-2 0-140 1-03 23-2 0-144 
: Table 2. Properties of materials used 0-069 9-20 10-6 0-132 0-97 24:8 0-150 
Material cand \Greund colt 0-069 6-53 85 0-106 0-782 17-6 0-113 
0-106 15-0 12-2 0-152 12 0-162 
ai Particle size . 90% bet- 95% bet 0-159 12-0 7-49 0-093 «60-686 14-0 0-092 
ween 06 ween 0-5 0-238 15-5 748 0-686 12-1 0-083 
0-238 11-9 6-13 0-076 0561 922 0-068 
Bulk density 0-156 19-4 788 860-098 00-500 13-6 0-130 
(Loosels | 10-0 5-66 0-071 0362 696 0-073 
packed) . . 1480 kg/m*| 560 kg/m’ | 0-238 31-9 12 14 0-137 l 
: Dynamic angle 0-238 5-0 6-02 O-O75 0-382 6-88 0-062 
of repose. . 36 40 0-156 20-0 7-63 0-095 0-226 40 0-103 
| 0-156 449 0-060 «0-143 2-57 0-058 
This table shows that the mater 0-238 582 0073 0-174 3-16 0-072 
ials were of rather uniform particle | (0-238 14-2 400 860-050 «0-119 209 0-048 
size. The density was taken for 
maximum percentage of voids which 0-046 1h2 990 0123 0-189 4-40 0-159 
0-077 4-75 0-059 0-091 2-26 0-080 
is probably nearest to the density 0-180 0-198 | 0-197 513 
actually occurring in the kiln. Since 0-180 2840 69.007 0-149 9.99 0-100 
both feed rate and hold-up were 0175 274 00384 0-:0526 0-029 
measured by weighing and ®, and X, (0-202 12-7 315 0039 00605 0-037 
are approximately proportional to 0-238 26-8 5-70 0-071 0-109 2-04 0-073 
each other (see Fig. 5), a small | 6 238 8-49 145 OO1I8 060279 0646 0-018 
error in density does not materially 
0-141 27-7 6-47 0-081 0-0632 1-30 0-085 
affect the experimental results. The | 0-140 10-2 2.40 0-080 00233 0-483 0-026 
angle of repose was measured during 0-0533 | 0-233 21-1 3.10 0-039 00302 0-602 0-033 
the experiments at different places | 0-238 12-9 1-81 0-023 O-0176 0360 0-017 
: in the kiln. No greater deviation 
0-048 10-6 6-80 0-058 00-0627 1-32 0-094 
than + 5% from the values indicated | 0-068 29-2 120 258239 
— | 26-0 750 0-004 0-069 1-45 0-106 
se the rotational speed of the kiln was 0-114 46-2 14-8 0-185 0-137 2.43 0-169 
sufficiently high to cause continuous 
movement of the material. Coke 
Table 3 summarizes the experi- 0-0100 0-083 2.98 2.8 0-092 =: 0-632 14-0 0-099 
mental results for the cylinder with- 0-0135 0-235 14-9 5 0-168 0-775 17-8 0-182 
out end constriction (h; = 0). As O0-O150 0-157 3-21 1-8 0-060 0-240 4-64 0-063 
© This piece of equipment was manu- 00232 0157 13-6 52 O172 0386 826 0-191 
factured by F. L. Smidth & Co., Copen- 0-0232 0-158 546 2-4 0-079 O177 3-28 0-081 
hagen, for the Laboratory of Chemical 0-033 0-235 3°57 090 860-0380 0-042 0-71 0-024 
Technology, Delft. It was kindly put at our 0-053 0-141 15-0 39 0129 0-095 193 0-145 


disposal by Professor H. 1. WATERMAN. 
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can be seen from Fig.5 both SuLtivan’s and our own 
results agree well with the theoretical curve, devia- 
tions greater than 30°% seldom occurring. In Fig.6 
the experimental and theoretical values of the total 
relative hold-up X, are compared. This figure clearly 
shows that there is some systematic deviation, the 
calculated value being somewhat too low for small X, 
and too high for great X,. In the region between 5 
and 10% hold-up, the average agreement is best. 
We also checked our results with SuLirvan’s 


empirical formula which can be written in the 


++ 


© SAND | mE ASUREMENTS 
a COKE 


x SULLIVAN'S DATA 
THEORY 


Fig. 5. General hold-up relationship for inclined kilns 


without end constriction. 


following form: 


2 
n tan « 


X (14) 


and where f is expressed in radians. For hold-up 
values greater than 5°, and in particular for small 
angles of inclination eq. (14) gives values for , a 
that are too high, sometimes by a factor 2. 

Two series of measurements were carried out with 
diaphragms at the discharge end of the cylinder, 
which had inner diameters of 0-177 and 0-158 m 
(h,/R = 0-10 and 0-20 respectively). The results 
are shown in Table 4. For purpose of comparison we 
calculated X, for the different experimental condi- 
tions, using first the theoretical line of Fig. 5 and 
correcting afterwards for h;, by means of Fig. 4. 
A comparison between the experimental and calcu- 
lated hold-up values is presented by Fig.7. This 
graph also includes the measurements of SULLIVAN 
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et al. with end constriction (tables 25, 26, 28 and 35 
of their paper). The agreement between (X,) 
and > 
the case h, 


cal 


is good, even somewhat better than for 
0 (Fig. 6). 


0,4 


0,3 t 


r/ exp 


oO 
+ + 
| \\ 
+ 
+ 
+ 
+ 


a COKE 
x SULLIVAN'S DATA 


+ 4 


0,01 0,02 0,03 004 0,06 0,08 0,1 0,2 03 04 


A. 


Fig. 6. Comparison between xperiumental and calculated 


hold-up values for kiln without end constriction. 


o SAND = 


a COKE 
x SULLIVAN'S DATA 
0,01 
0,01 0.02 0,03 0,04 0.06 0,08 01 0.2 0.3 04 


Fig. 7. Comparison between experimental and calculated 


hold-up values for kilns with end constrictions. 


4. Discussion or THE THEORY 
It is rather surprising that such a simple model of 
operation as described in section 2 yields results 
which appear to be useful for practical design work. 


However, the actual motion of the solids is much 
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Table 4. Experimental results for hy R 


0-10 and 0-20 
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Engineering Science 


throughout this process. We believe 


Hold-up ht 
7) No 


Sand, hy, R 


10-8 O-OL7 


6-0 O-058 
11-7 
3-0 


0-0230 57 aS 


O-OL35S 0-022 


00-0230 0-027 
O-O78 
O-OD35 O-OS 1 


San 


0-0100 0-153 0-037 


0-156 0-102 0-108 


Morse O-238 0-140 O47 


0-134 


0-36 


0-0240 


508 


O-156 


O-156 O-O87 


O26 


0-141 o-102 


O-OL00 O-TS5 


0-127 0-026 


O-OL35 0-016 


00-0230 0-204 


O-O76 


0-038 


O-O230 230 
1-45 


4-00) 


0-0535 0-132 
hy 020 
3-10 
2-35 
4-58 


00-0230 2-4 aS 


OOTT 
O-151 
0-191 


O78 
0-152 
056 
O4M6 
0-0230 2.00 0-20 
1-30 


$45 


O-O4S “36 


0525 0-147 OOT74 


that the assumption of a flat sur- 
Vo N 
—s face introduces no serious error, as 


long as the average angle of 


repose under actual conditions of 


0-160 operation is used, and the centri- 
0-069 
0-132 
0-164 


fugal acceleration caused by rota- 


tion is negligible with respect to 


gravity. 

The assumption of an infinite 
O-065 
cascading velocity of the particles 
might, on the other hand, seriously 
0-156 affect the results. A limited average 
velocity causes the calculated hold- 
the 


0-159 up to be too small, because 


mos 


laver of material flowing down re- 
presents an additional hold-up not 
taken by the theory. In 


em 


care of 
separate experiments in a 


diameter cylinder, we measured 

superficial velocities of the order of 
0-191 
140 


0-199 


In this case the material 
this 


about 20° of the total hold-up. 


30 cm sec. 


involved in movement was 

i From the mechaniom of flow, which 
O73 
0-044 


0-154 


is not unlike that of liquid films, 
this percentage may be expected 
kiln 


ascertained 


with increasing 
This 
measurements in 


evlinders of 2-4, 46 


to clecrease 
diameter. was 


0-102 by three small 


0-070 and 11-7 em 
0-139 


diameter at constant A/R values. 


O58 
0-033 


0-178 


b) Gradient along the kiln 


From eq. (2) which has been derived 


for small gradients, it follows that 


more complicated than would appear from the 


assumptions made. Moreover, a number of simpli- 
fications have been introduced in the calculations, the 
significance of which for the final results cannot easily 


he seen. A few of these points will be discussed here. 


a) Movement of the particles 

In reality, the surface of a granular material in a rotat- 
ing cylinder is not flat but wave-like. This form of 
surface is connected with the fact that, once a 
particle has reached the surface it first is accelerated, 


then reaches a more or less constant velocity and 


finally slows down near the bottom of the slope. 
that 


It appears the angle of “repose’’ decreases 


dh/dx can become very great for 
relatively small values of 4/R. In reality this value 


cannot exceed tan For cases where hp = 0 and 
the end effect is relatively important (N, N,> 1) this 
will lead to calculated hold-up values which are too 


high. 

On the other hand, for kilns with a thin diaphragm 
at the discharge side, as used in our experiments, the 
calculated hold-up will be too low. The constriction 
gives rise to a very complicated flow pattern. It 
results in a positive contribution to the theoretical 

dh/dx value, due to the restircted passage. The 
latter effect the 


becomes less important. 


also decreases as end correction 
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c) Simplifications used 


Whereas the approximation of eq. (8) is sufficiently 


accurate, the substitution eq. (3) is rather crude, 
especially for low values of 4/R. It would cause the 
lowest calculated hold-up values to be too high. 
However, this approximation has the important 
effect of reducing the number of variables in eq. (2) 
from 3 to 2, which made a rather simple representation 
possible. Then, however, a close correlation between 
the two variables of Fig.5 can no longer be expected. 
Whereas the experimental results are estimated to 
have an accuracy within 5%, the mean square 
error of all the + ae (X,).4, Values amounts to 
about 

It appeared to us impossible to take into account 
any of the effects mentioned without making the 
calculations extremely complicated. The experimental 
results indicate that the effects of the various as 
sumptions and simplifications have, to a great extent, 
cancelled each other out, leaving only a slight SVs- 
tematic deviation as shown in Figs. 6 and 7. Whether. 
however, these fortunate conditions persist also for 
industrial kiln dimensions has to be seen from ex- 
perience. The following example, taken from SAEMAN’s 
paper gives an indication that the theoretical results 
are applicable to large scale equipment. 

Example. An industrial nodulizing kiln consists of two 
sections, Section | on the discharge side has a length 
L $5 ft (13-7 m) and radius RK, ft (005m). Of 


1 
4 SO ft (24-4 m und tt 


section 2 on the feed end, L, 
(107 m). Further data are 
— 6&1 cu ft min (2-88 1073 sec) 
45 


rev sec 


tan x 5 in ft (0.042) 
Measured time of passaye T 40 min. 
From these data we find for section 1: 
Ne, 115 
Ney 1-17 | 


According to Fig. 5 this value corresponds to X Nv; ‘= (O70, 


No, 1-35. 


giving \, — 0-0565 and the time of passage in section |: 
12-6 min. 
In order to calculate the hold-up of section 2, the height of 
the bed A, at the junction of the two secticns has first to 
L 
be calculated. Now 
Ay R, R, +h, 


where A, is the height of material at the entrance of section 1. 
From Fig. 2 we find h, R, Ny, 


0-190, giving = 0-21 m 


and hy 0-33 m. Thus for section 2: 
hy No O-SO5 


} O-387 


No, Nj; 0-60. 
/2 Nx, 0-742 | 
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965 


3 


Using these values and Fig. 4 and 5, we find X, N; 0-035, 
X,. = 0-052 and so for the time of passage in section 2: 


tT, = 26-2 min. 
The total time of passage 1, +1, 
with the 40 min actually observed. 
The authors wish to thank Messrs. L. E. M. 
Muppe, E. F. Bunge, A. Strxker and A. E. F. 


CAPPELLE, who assisted in parts of this investigation. 


38-8 min agrees closely 


NOTATION 
All quantities are expressed in units of the practical system 
(m, kg-mass, sec). 
f; correction factor for end constriction with 
height h, 
h depth of bed of material (see Fig. 1 
OW hy item at feed end and discharge end respectively 
hk length of path of cascading particle perpendicula: 
to kiln axis (see Fig. | 


L length of the kiln 


n rotary speed of the kiln 

R radius of circular cross-section of the kiln 
$ axial component of the path of cascading particle 
n axial distance from feed end 

X fraction of cross-section taken in by material, 


dimensionless 


X, = average of X, over kiln length, relative hold-up, 
dimensionlk Ss 

% angle between kiln axis and horizontal plane 
jh angle of repose of material 

4 angle on path & (see Fig. 1) 
Pp, volumetric flow of material 

T time of passage through kiln 
Np sin p n tan «, dimensionless 

R cos L tan x, dimensionless. 


Note Appep IN Proot 
Recently, our attention was drawn to a paper Picke- 
RING, Feakes and Frrzceravp, “Time for passage of 
material through rotary kilns”’ (J. Appl. Chem. 1951 
1 13). The adopted transport mechanism is the same 
as in the present paper but the influence of the varying 
height of the material along the kiln has been neglected. 
Whereas these authors find a systematic discrepancy 
of about 10° between their simplified formula and 
their experimental results, the design equations 
given in the present paper reduce this difference to 


2"> on the average. 
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a A new method for predicting the plate efficiency of bubble-cap columns 
Le 
S. BAKOWSKI 
‘ag Imperial Chemical Industries Limited, General Chemicals Division, Research Department, Widnes 
Received Auqust 1952) 
Summary Recent investigation of the mechanism of bubble formation at submerged slots has shown that 
“ under the conditions normally existing distillation columns the vapour may be regarded as flowing 
a through channels, rather than passing im the form of streams of separate bubbles 
’ By a simplified treatment based o1 mechanism an equation has been derived for calculating the 
a local efficienev of a bubbk cap |} late | effimencies calculated by means of this equation have been 
ompared with published experimental data for columns ranging from small pilot-plant units to large stills 
| 
of several feet diameter as used in the o lustry; data for absorption processes in which the gas film offers 


the principal resistance to 
the predicts d and measured efficiencies 
rhe difference between the local effi vy and the 
usually to be small, and the new equatior 
A sper 


indicated 


Résumé— Des recherches récentes sur 
que dans ! 
de canaux continus plutot que sous fort i 
sappuyant sur cette hypothes 


locale d'un plateu A calottes Les valeurs 


diverses colonnes (colonnes de petites in 


plusieurs pieds de diamétre). L’ auteur ex 
ou le film gazeux constitue la résistance pri 
entre les efficacités calculées et celles qu sol 
Pour des plateaux réels, il y a peu de différen 
la nouvelle équation peut étre utilisée en pr 
L’auteur envisage une application particuliére de 


distillation dans le vide 


matter transt have also been examined 


be used im practice 


ial application of the new equation to the 


tormation cle 


s conditions normales des colonnes de distillation 


ipale au transfert de matiére 


Good agreement was found between 


Mcurruree plate-efficiency for actual plates appears 
to predict the MurPHREE plate efficiency 


design of plate columns for vacuum distillation is 


bulles au travers de fentes immergées ont montr 


écoulement de la vapeur s'effe ctue sous forme 


hapelets de bulles séparées 
vuteur a établi une équation permettant le caleul de lefficacit 
ainsi obtenues ont été comparées aux résultats publiés pour 
tallations pilotes jusqu’aux colonnes de l'industrie du petrole ce 


également les données relatives A des operations absorption 


Il vy a un accord satisfaisant 


t mesurees 
e entre lefficacité locale et l'efficacité suivant MURPHREE; 
itique pour prédire lefficacité de 


son équation au caleul des colonnes & plateaux pour 


One of the main problems in the design of a_ efficiency and the product of the relative volatility 
;. bubble-plate column is the proper estimation of and viscosity of the mixture distilled. According to 
* plate efficiency. Three different meanings have been O’CoNNELL the plate efficiency decreases with increase 
attached to this term, when defined in terms of the in the value of this product. This method gives 
‘ Vapour composition the most commonly used over- fairly good agreement in some cases when applied to 
. all efficiency” is the ratio between the number of distillation of hydrocarbons, but, for instance, the 
. theoretical plates and the actual number of plates efficiencies predicted for aleohol/water mixtures were 
required to effect a given separation of the compo- found to be considerably lower than those encountered 
nents; the McurPHREE efficiency’ is the ratio between in practice. 

the actual enrichment of vapour produced by a plate A method based on consideration of the diffusion 
and that which would be obtained if equilibrium occurring between individual vapour bubbles and the 
between liquid and vapour were reached. If referred liquid over the slots was proposed by Geppes [9] 
to a complete plate in a column the latter is termed and recently developed by Ju Cars Cuv, Donovan, 
the “Murpuree plate efficiency’’, and if referred to Bosewettand Furmersrer [5]. Although this method 
Z a single slot or point on the plate it is termed the gives reasonably good agreement with several ex- 
: “local” or “point” efficiency. Lewis [20] derived perimental results reported in the literature, critical 
“s analytical expressions showing the relationship be- examination shows that the assumptions made by 
tween the plate and local efficiencies for three theoretic- | Geppes diverge considerably from the actual condi- 
7 ally possible cases tions in a distillation column. This becomes clear if 
A great deal of experimental work has been done comparison is made between the high vapour rates 
and several methods for predicting plate efficiencies on the plates and the small size and low velocities of 

; have been proposed. One of these, originally proposed the bubbles assumed in this treatment. 
by DrickaMerR and Braprorp [7] and subsequently The mechanism of bubble formation at submerged 
modified and developed by O’ConNneLL [25], is based slots was recently investigated by Sretis and Ba- 
: on an empirical relationship between the overall Kowsk1 [30]. A new method of predicting plate 
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efficiencies based on the results of these experiments 
forms the subject of this paper. 


MECHANISM OF Vapour FLOow LiguIp on 

A PLATE 
In preliminary experiments [30] in which bubble- 
formation at single and multiple slots was studied 
with the aid of stroboscopic illumination it was found 
that the properties of the gas and liquid and the shape 
of the slots had little effect on the apparent size of 
the bubbles; the main factors were found to be the 
rate of flow and the slot submergence. There was 
frequently a discrepancy between the known gas 
flow-rate and that calculated trom the size and fre- 
quency of the bubbles, particularly tor multiple slots, 
the former being the greater: this effect was most 
prevalent for multiple slots for which the bubbles were 
smaller. This suggested that intermittent channels 
were formed, connecting the slots directly to the 
surface of the liquid. 

Sretts and Bakowsk1 [30] investigated the 
mechanism further for single slots using high-speed 
cinematography. It was confirmed that at flow rates 
and slot submergences normally employed in distilla- 
tion columns there are few intervals between individual 
bubbles, these being connected with each other more 
or less continuously so that direct flow of gas from 
the slot to the surface of the liquid can take place. 
Further work by the same authors, not yet published, 
has revealed the existence of similar continuous 
channels when gas bubbles are formed at closely- 
spaced multiple slots as used in distillation columns. 

It is clear that under these conditions there is no 
simple relation between the volume of an individual 
bubble, the rate of flow, and the time of bubble 
formation. Each bubble can be regarded as a channel 
connecting a slot with the surface of the liquid. 

The rate of vapour flow through each slot on 
a bubble plate is directly proportional to the super- 
ficial vapour velocity in the column, and therefore 
to the boil-up rate or vapour load. It is usual to 
maintain this at the maximum value allowable without 
causing excessive entrainment of liquid, in order to 
reduce the diameter of the column as far as possible. 
Safe values of the vapour load, which depend on the 
densities of the liquid and vapour and on the plate 
spacing and depth of liquid seal, can be calculated by 
one of the conventional methods [2]. The correspond- 
ing vapour rate per slot can then be derived for any 


design of bubble plate. It is evident from Fig. 1 which 


A new method for predicting the plate efficiency of bubble-cap columns 


shows diagrammatically a common arrangement of 
3-in. diameter bubble caps at 4}-in. triangular pitch, 
each having 38 rectangular slots, that the horizontal 
expansion of the bubbles is limited by those issuing 
from neighbouring slots. Thus in the case illustrated 
the free gas passage would be limited to about 0-9¢em?* 
and for a superficial vapour velocity of 60 em sec 
in the column the linear velocity through this cross 


section would be of the order of 200 cm/sec. This 


Fig. 1. Typical arrangement of bubble caps. Lateral 


growth of the bubbles is limited by the bubbles issuing 
from neighbouring slots. 

far exceeds the normal velocities of gas bubbles rising 
freely in a liquid, and thus also suggests that channels 
must be formed. 

Although increasing the slot submergence has 
a favourable effect on plate efficiency this should not 
be carried too far because of the consequent increase 
in liquid hold-up and pressure drop. In practice the 
total slot 
5em, and a }-in, (1-27 cm) static seal is often re- 


submergence does not usually exceed 
commended as normal for distillation at atmospheric 
pressure [6]. Under these conditions, i.e. high rates 
of flow and comparatively low slot-submergence, the 
formation of separate bubbles is improbable, and the 
most reasonable assumption seems to be that the 
vapour flows through channels which may be re- 
garded as short wetted-wall columns in parallel. It 
will be seen from Fig. 1 that area of the free liquid 
surface may be expected to be small in comparison 
with that presented by the walls of the channels. 
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PREDICTION OF THE LocAL ErriciENcy BASED ON THE 
CHANNELLING MECHANISM 


The process of channelling of vapour through the 


somewhat complicated and the 


A simplified 


liquid is actually 
channels are of irregular cross-section 
mechanism in which the channels are regarded as 
wetted-wall columns of uniform cross-section was there- 
fore assumed as a basis for further treatment, as 
illustrated in Fig. 2 

The treatment involves the following assumptions 

a The channel length Is equal to the total sub 


mergence, 1.¢. to the sum of the liquid seal over the 


Tike alone al ‘ hanme lling hanian 


Fig » 


slots and the active slot opening, given bv the « lation 


The mean cross-sectional area of the « snnels 
regarded as approximately constan er the 
slot 

annel 


The rate of matter transfer through the 


walls is governed by the resistance of the film 
cd) The effects of temperature and of the properties 

of the Vapour on the matter-transfer coefficie 

be expressed by a convenient approximation 

namely tl 


on the elementary kinetic theory 


K, x (MT) 


where WV is the molecular weight 
T is the absolute temperature 
which concern the 


dimensions of the 


first two assumptions 


channels are based on actual 
measurements of the volumes of bubbles formed at 


slots 


] De rivation of the qe ne ral equation 


The local or point efficiency for a binary mixture in 


terms of the partial pressure of the more volatile 


licting the plate efficiency of bubble cap columns 
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component, as defined by Murpures, is given by 


the equation 
(2) 
where partial pressure of more-volatile com po- 
nent In vapour entering the channel 
(from plate below) 
partial pressure of more-volatile com- 
ponent in vapour leaving the channel 
partial pressure of more-volatile com. 
ponent in equilibrium with the liquid, 
Considering the vapour issuing from a single slot: 
since approximately equimolar quantities of the com. 
ponents are exchanged between the liquid and vapour 
the vapour rate may be taken as constant and the 
balance can be 


matter-transfer expressed by the 


equation 


273 
/ K 


22400) 


A Ip (33) 


m 


where vapour rate per slot (em®* sec) 
absolute temperature (°K) 
matter-transfer coefficient 

(g mol sec . em? atm) 

log mean driving force (atm) 
surface area of channel (cm*) 


yviven by 


and substituting this in eq. (3) gives 
AT 


In (5) 


Combining this with eq. (2 then leads to the following 
equation tor the local efficiency 
AL AT 


9731 (9) 


In order to apply this equation in practice it is 
necessary to express the matter-transfer coefficient 
and surface area in terms of the operating conditions 
and the properties of the substances involved. 

Taking the rate of matter transfer to be controlled 
by the gas film the coefficient may be considered to 
vary with the vapour velocity and the pressure accord. 
ing to an equation similar to that given by Carey 
for a wetted-wall column 

P 


ky pow gon (4) 


and [4 


For equimolar counter-diffusion, however, the drift 
term P p, is unity and can be omitted, so that: 


K ky pees (8) 


J 
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where k, is a constant depending on the temperature 
and the properties of the vapour 
v is the linear velocity of vapour in the channel 
(cm sec) 
d is the mean equivalent diameter of the 
channel (cm) 
P is the total pressure (atm). 

If the rate of matter transfer is assumed to be 
proportional to the difference between the numbers of 
molecules colliding with unit surface in unit time when 
under the actual and equilibrium partial pressures, it 
can be shown from the elementary kinetic theory that 

273 
22400 Mee 
where M and 7 


temperature (“K). This relation has been successfully 


(9) 


molecular weight absolute 
used for calculating rates of matter transfer in high- 
vacuum evaporation and molecular distillation, and 
may be expected to hold when the yas film thickness 
and mean free path of vapour molecules are of the 
same order of magnitude. It is in reasonable agree 
ment with the experimental data of JounsTone and 
PicrorD [15] in a wetted-wall column. It is conveni- 
ent to employ this simple relationship in preference 
to a complex function of viscosity, density and dif- 
fusion coefficient as used in the normal treatment ot 
matter transfer, using suitable mean values for the 
molecular weight. 

Introducing now the assumptions based on the sim 
plified channelling mechanism, namely that the channel! 
length A is equal toh hand that the mean cross-sec- 
tional area of the channel is constant it follows that 

Vapour velox the rate per slot 

Interfacial area A x h. 
Equivalent diameter — constant. 

Combining eqs. (7), (8) and (9), taking d as constant 
and substituting values of » and A in terms of d there 
is obtained on conversion to common logarithms: 

l kh 7 


where & is a constant to be evaluated by comparison 


with experimental data 


” Evaluation of th empirical constant 


Preliminary examination of data for a number of 
distillation processes showed that reasonably constant 
values of k were obtained. The data of Carry, Gris- 
woLpD, Lewis and McApams [3] for the systems etha- 
15 experiments 


and benzene toluene in 


nol water 
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with rectangular slots of normal width } in. (0-32cm) 
and }in. (0-64 em) gave a mean value of k = 0-34. 
The final equation for the local efficiency therefore 

becomes: 
l 0-344 7° 


log io 
COMPARISON WITH EXPERIMENTAL DATA 

Published experimental data on plate efficiencies for 
a range of systems and types of column have been 
compared with the efficiencies predicted by eq. (11). 
The constructional details of the columns, the operat - 
ing conditions and the systems employed are given 
together with the observed and predicted efficien- 
Table 1. The 


trated graphically in Figs. 3 to 6, 


cies in agreement obtained is illus- 


which show the 
expermental figures together with lines calculated 
from eq. { 11) 

The methods used in evaluating the variables in 
eq. (11) for the experimental conditions and in making 
the comparison are described in the following sections 


of the paper, 


I. Vapour rate per slot, V, 


The flow rate per slot was computed by dividing the total 


vapour rate in cm®/sec, calculated from the superficial vapour 


velocity and the area of cross-section of the column, by the 


number of slots per plate. 


= Mean molecular weight, M 


This was calculated from the composition of the mixture in- 


volved. Since at total reflux the cor position of the vapour 
entering the plate is equal t that of the liquid on the plate 
the latter composition was used for computing WM. When the 
concentrations of the components were not given MW was taken 


as the arithmetic mean of those of the components. 


é. Total slot subme rence h 


The total submergence / was obtained by adding to the static 
calculated by a modi- 
This 

the 


liquid seal h, the active slot opening A, 


fication of the equations of Rogers und Taree [28]. 


modification consisted of neglecting the density of 


vapour 9, compared with that of the liquid (i.¢. assuming 


0 0 o,) and expressing the vapour density as a function 


of pressure, temperature, and molecular weight. The resulting 
equations are as follows: 


For a rectangular slot: 


h, = 0-038 (M o,)' (V,/W) (12) 
For a triangular slot: 
h, = 0-2 (M P/To,)° (H VW) (13) 


where VW width of slot at the base 
H height of slot. 
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For trapezoidal slots it was necessary to solve graphically the 
equation: 
To, \° (W 


DD PM 2-5 h' 


Wp) 


H (14) 


where W, and W are the widths at the top and bottom of the 
slot respectively. 

At high rates of flow it was sometimes found that the 
calculated active slot opening exceeded the total height of the 
if the 


caps rested directly on the plates 4, was taken to be the slot 


slot. In such cases, slots were closed by a skirt or the 

height 1/7; otherwise 4, was corrected for flow of vapour beneath 

the cap. 
Although the 


head of liquid over the exit weir Peavy and Baker [26] have 


slot submergence may be increased by the 


shown that liquid is often pumped over by the vapour-lift 
action of bubbles from neighbouring caps and its level can 
even be somewhat below the weir. In general therefore only 
the static liquid seal was taken into account in computing the 
total submergence 

In sore however. a corrected efficiency E’ was also 


ileulated using as slot submergence 


where 4, is the liquid head over the weir, « lated from the 


well-known Francis formula. 


/ (‘om parison ol obse rv d and pre dict d CTTICIOHCIES 


Since the value of the constant in eq. (11) was based on ex 
periments in a small column with a single cap on each plate 
In large columns 
higher 


mixing. In 


it should give the local or point efficiency. 


the Murruree plate efficiency should be owing to 


cross-flow with incomplet some 


CASES 


licguid 
where efficiencies greater than unity were observed for instance 
when recalculating the data of Peavy and Baker [26], an 
ittempt was made to derive the local efficiency from the ex 
perimental data by means of Lewis’ eq. (20) for Case I 
(assuming complete mixing of vapour and no mixing of liquid), 


EyyvimV/L) 


In 


where E, 
m — slope of equilibrium line (dy*/d x) 


V/L 


local efficiency 


Murruree plate efficiency 


ratio of molar vapour and liquid rates. 

In general, however, the differences between the observed 
values of Ey, and the point efficiencies calculated from 
eq. (11) were found to be small and much lower than would 


be expected from Lewis’ equations, 


DiscuUSSION 
It will be seen from Table 1 and Figs. 3 to 6 that in 
most cases the efficiencies calculated from eq. (11) 


agree well with those determined experimentally. 


A new method for predicting the plate efficiency of bubble-cap columns 


From the results it has been found possible to evaluate 
the relative importance of the main factors which 
affect the efficiency of bubble-cap plates: these are 


the slot dimensions; the properties and composition 


VAPOUR RATE sec SLOT 


Variation of efficiency with vapour rate, 


width and liquid seal. 


ALCOHOL-WATEF 


FERICIENCY 


Fig. 4. Variation of efficiency with vapour rate and liquid 


seal. 


of the mixture handled 
the 


the operating pressure and 


temperature ; vapour rate, and the slot sub- 


mergence. 


1. Effect of slot width 


The width of the slots seems to have no direct influence 


on plate efficiency; but there is an indirect effect 


caused by the change in active slot opening and total 


| 
2 
wore cm 
0% 07% 063 26 
in 
ao} 
| 
70 } 
7 
| 
| | 
=THYL ALCOHOL WATER | 
an 
} 
| 
OL. } 
‘ 
‘ d - 
52 +h +h, 
~ hy = 18cm 
70 
fo) 
are fo} 
a 
AC 
a 
k 1} 
J 
- 
0 400 200 300 400 
Vapour Rare, Cm’/Sec PER SLOT 
-- 
‘6 
277 
3 
‘ 


Chemical 


8S. Baxowski: A new method for predicting the plate efficiency of bubble-cap columns Engineering Science 
or 
| FQ 4 
0} 
| 
°° 
° 
| 
| © 
| ¢ 
(ZENE Benzene - To.vene 
| 
u 
| hy = 1-1 
©| 
300 400 fon ( 400 200 300 
apoun Rare. cm/sec Per VAPOUR Rat PER SLOT 
Fig. 5a Fig. 5b 


EFFICIENCY 
8 


© PAGE § WICKER’S DATA hy =_0.7 om 
ae 
CAREY'S DATA 
L ] 
00 200 300 400 500 50 AC 200 
VAPOUR RATE, Cm5/Sec PER SLOT VAPOUR RATE om? Sec PER SLOT 
I Fig rd 
an Zz 
Ea 
+ 


o 


& 

5 

iS 


o/ 
i] 


> 
30 + + 
2 
“a: 2 
wi 


- 200 400 600 800 900 
Vapour Rare Cm’/ sec PeR Sor VAPOUR RATE Cm3/Ser Der 


Fig. 5e. Fig. 5f. 
Figs. 5a—f. Variation of efficiency with vapour rate. 
278 


= 
4 
90 
A | 
x ac VOL. 
E* 
| 
x 
|_| 
BENZENE - TOLUENE = 0.30m TRICHLOROE THY LENE | 
a 
j 
— WATER 0} ALCUMOL —— WATER, 


8. Bakowskr: A new method for predicting the plate efficiency of bubble-cap columns 


] 
| 
| 


LUENE 


6 
ACTION BENZENE 


Fig. 6a. 


WATER 


- 3:0 


EFFICIENCY 


t—J 
|__| 
| 
| 
| 
| 


AZeorropic Mixture 


4 


| 


3 4 ¢ 
Fraction 
Fig. 6c. 
Figs. 6a—c. Variation of efficiency with vapour 


composition 


submergence at a given vapour rate. This is evident 
from the work of Cargy ef al. (3) (Fig. 3) in which 
The 


pitch of the slots in these experiments varied from 


the slot width was varied from 1-6 to 12-5 mm. 


6 to 48 mm: even under these conditions the assump- 
tion that the vapour passes in the form of channels 
rather than as separate bubbles appears to be justified. 


2. Effect of properties and composition of the mixture 
The effect of the properties of the mixture was studied 
by Gapwa [8]; except for mixtures of methanol and 
isobutanol his results agree fairly well with the pre- 
dictions of eq. (11). More recently experiments with 


various liquid mixtures were carried out by WILLIAMs, 


EXCERWENTAL CuRve 


BENZENE 


TOLUENE 


5 6 


MOL. FRACTION BENZENE 


Fig. 6b. 


STIGGER and NicHois [32], who, however, reported 
the results as overall efficiencies. These results are 
unfortunately too scattered for a comparison with 
eq. (11) to be made, possibly because the liquid seal 
employed in this work was low. 

Eq. (11) implies that molecular weight and abso- 
lute temperature affect the plate efficiency. Since 
liquid viscosity decreases as the temperature rises and 
usually increases with molecular weight this may ex- 
plain the relationship between plate efficiency and 
viscosity which forms the basis of the correlations of 
DRICKAMER and Braprorp [7] and O’ConneELL [25] 
for overall efficiency. 

The effect of the composition of the mixture has 
been studied by Carey and Griswoip [3] for the 
henzene toluene system (see Figs. 6a and b); it is 
evident that the effect is small, but it would be ex- 
pected to be greater when the molecular weights of 
the components differ widely 

The experiments of LanGpon and Keyes [19| 
with isopropanol/water show good agreement with 
eq. (11) except when the composition is close to that 
of the azeotrope (see Fig. 6c). Similar low efficiencies 
were also reported by Keyes and Bymawn [17] for 
ethanol water when close to the azeotropic composition 
and also for very weak alcohol mixtures. These 
discrepancies may be attributed to lack of precision 
in determining the Murpuree efficiency in these 
ranges of composition, the results there being very 
sensitive to experimental errors. 

Examination of data for air humidification and 
absorption of ammonia from air into water [10], [31] 
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has also shown that eq. (11) gives reasonably good 
results when applied to absorption and desorption 
processes in which the gas-film resistance is the limiting 
factor. 
3. Effect of pressure and temperature 
The effect of pressure on plate efficiency can be seen 
in Table 1 for the experiments on air humidification 
31] and for tests on columns used in the petroleum 
industry. Vacuum distillation experiments by Kirscu- 
BAUM [18] on the system ethanol water also indicate 
increased efficiencies at the lower pressures 

At a given vapour rate per slot (em® sec) the direct 
negative effect of increased pressure 1s COM] nsated 
to some extent by the increase in active slot opening 
h,) with pressure [eq. (12) and (13) 

The predicted variation with temperature, which 
appears m eq. (11) as the ratio (T/M)**, cannot be 
verified from the data examined since the range of 


values of this ratio covered is too small 


4. Effect of vapour rate 

The effect of vapour rate has been studied by several 
observers in order to establish the optimum operating 
conditions. At vapour rates above some crit! il value 
depending on the plate spacing entrainment of liquid 
becomes sufficient to reduce the efficiency. The gene- 
ral conclusion which can be drawn from the published 
data is that except at very low rates (uneconomic In 
practice) or at vapour velocities high enough to cause 
excessive entrainment the plate efficiency varies little 
with vapour rate. 


This conclusion is consistent with eq. (11 which 


0-25 


that log |1 (1 


so that a small decrease in efficiency with vapour rate 


shows E)) is proportional to 
per slot should occur. This direct effect is « ompensated 
by the increase in active slot opening, which gives an 
increase in total slot submergence. When the static 
seal (h,) is low the latter effect should prevail, so that 
under these conditions there may even be a slight 
increase in efficiency with vapour rate A rapid drop 


in efficiency would be expected for vapour rates 
exceeding that at which the whole slot opening just 


becomes active 


>. Effect of slot submergence 

The most important factor influencing plate efficiency 
is the slot submergence, which depends on the static 
seal prov ided. At very low submergences quite small 
variations in liquid level may considerably affect the 


efficiency; under such conditions the experimental 
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data are somewhat scattered, as may be seen in the 
Figures 5b and e. 

In large columns differences in liquid seal resulting 
from hydraulic gradient or from faulty levelling of 
the plates may thus directly affect the efficiency, apart 
from the reduction caused by the accompanying 
maldistribution of vapour between unequally sub 
merged caps 

The method outlined in this paper for predicting 
plate efficiency involves a relation between efficiency 
and slot submergence, and is therefore considered to 
be superior to that of O’CoNNELL, which does not 
take account ot this factor. A modification of O’Con 
NELL’s method incorporating an allowance for liquid 
seal has however been made recently by Ju Carn Crt 


et al 1D 


APPLICATION OF THE New Equation To THE DeEsiGn 
or CoLUMNS 
Theoretically eq. (11) should give the local efficiency 
and the Murpuresr efficiency for whole plates should be 
somewhat higher. As the observed MuRPHREE plate effi 
ciencies differ little from the values predicted ; eq. (11) 
should be used in practice to predict the MURPHREE 

plate efficiency to be expected from a column. 
This is not necessarily equal to the overall efficiency 
defined as the ratio of the number of theoretical plates 
to the actual number of plates required. They are 
equal only when the equilibrium and operating lines 
in a McCapne-Tutece diagram are parallel, i.e. when 
l. 


to employ the Murpuree efficiency in conjunction 


In important cases of design it is preferable 


with the McCane-TureLe diagram to obtain directly 
the number of actual plates necessary. This may be 
done as illustrated in Fig. 7 by drawing the vertical 
steps so that the enrichment per plate is Ey; times 
that corresponding to equilibrium. This is most 
easily done by starting from the reboiler end of the 
column, but this is not essential. 

A special application of the new method may be 
made in the design of columns for vacuum distillation 
of heat-sensitive substances, in which the pressure-drop 
is often limited by the temperature permissible in the 
reboiler. The pressure drop per plate can be lowered 
by reducing the liquid seal, but this in turn reduces 
the efficiency so that more plates are required. If 
a mean value of the overall efficiency equal to » times 
Eyy, can be taken, the new method can be used to 
determine the best combination of liquid seal and 


number of plates, as follows: 
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Let ., — the number of theoretical plates required 
in the column 

N <= number of actual plates 

Z — specified pressure drop expressed in em 

of liquid. 
Then, neglecting pressure drop in the caps and 
risers 
Ny Z 


nEuv M+he (15) 


But by eq. (11), taking Fy), E,, and P as the mean 
absolute pressure : 
logis l 0-34 (Ar (11) 
l Ey Mes pes Vow 
A value of A, can be found by trial such that both 


these equations give the same value for the efficiency : 


E A. 


VAPOUR 


FRACTION OF Low BOILER tt 


Ma 


Xw Xp 
MOL FRACTION OF Low BOILER IN LiquD 


Fig. 7. Modification of diagram to allow 
for Murruree plate efficiency. The vertical steps are 


IE 
taken so that AS wi 
At Dk 


insertion of this value in eq. (15) then gives the number 
of actual plates. If a negative value of h, is found this 
indicates that the specified pressure drop is too low 
and the pressure at the column head will have to be 


reduced if the desired separation is to obtained. 


CONCLUSIONS 

1. Recent investigations of the mechanism of bubble 
formation at submerged slots have shown that under 
the conditions normally existing in distillation columns 
the vapour may be regarded as flowing through a chan- 
nel, rather than passing in the form of a stream of 
separate bubbles. 

2. By a simplified treatment based on this mech- 


anism an equation has been derived for predicting 
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the plate efficiency in bubble-cap columns. The main 
factor affecting the efficiency is slot submergence: 
other important factors are vapour rate, pressure, 
temperature, and the molecular weights of the com- 
ponents. 

3. Comparison of experimental data published 
in the literature with efficiencies calculated from the 
new equation shows good agrement, apart from a few 
discrepancies at the extremes at the concentration 
range and in the neighbourhood of azeotropic mixtures. 
These are probably due to the lack of precision inherent 
in the determination of the MuRPHREE efficiency in 
these regions of composition. 

4. The method can also be used for absorption 
processes in which the gas-film controls the rate of 


absorption. 
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NOTATION 


1 interfacial area for matter transfer (cm*) 
d — equivalent diameter of channel (cm) 
E yy) MvurpPureE plate efficiency (based on the change 


of vapour composition) 


E,, = overall efficiency (the ratio of the numbers of 
theoretical and actual plates required) 
local or point efficien 


H Total slot height (cm) 
h total slot submergence (hy h.) (em) 
total submergence corrected for weir head (A,,.) (em) 


h, — static liquid seal (cm) 


h, — active slot opening (cm) 
h,, head of liquid over weir (cm) 
kK overall matter transfer coefficient 


g mol/sec . cm* atm) 


gas-film coefficient (g mol/sec cm*® atm) 

L Molar liquid flow (g mol/sec) 
M molecular weight 

m — slope of equilibrium line (d y*/dx) 

N number of theoretical plates 

\, = number of actual plates 

P — absolute pressure (atm) 

p — partial pressure of more volatile component (atm) 


T absolute temperature (“K) 


1’ = molar vapour flow (g mol sec) 
\. = vapour rate per slot (cm* sec) 
v = vapour velocity in channel (cm sec) 


WW slot width (cm) 
W, = width at top of trapezoidal slot (em) 
r mol. fraction in liquid 


y — mol, fraction in vapour 


281 


‘ 
, 

| 7 
| 
952 — 
| 
| { 
— 


Chemical 
E. Wieke and W. Votitmer: Stromung von Gasen durch feinporige Stoffe Engineering Scien 


y* equilibrium mol. fraction in vapour Inst. Chem. Eng. 1946 41 393 13| Grisworp, J. and ST 
wart, P.B.; Ind. Eng. Chem. 1947 39 753 14) GUNNESS 
R. C.; Ind. Eng. Chem. 1937 29 15!) Jounstrone, H. 


and Prororp, R. L.: Trans. Amer. Inst. Chem. Eng. 1942 38 
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Strémung von Gasen durch feinporige Stoffe 
(Flow of gases through micropores) 
E. Wieke and W. Vo_umer 


Physico-chemical Institute of Gottingen Universit, 


Summary The flow of gases in the transition zone between laminar (PorseutLLe) and molecular (KNUDSEN 
flow is of great importance in vacuum technique is finding increasing application in modern chemical 
industry. Moreover, there are special separation effects in this transition zone*. Our knowledge of the 
flow through micropores, which often takes place even at atmospheric pressure in this transition zone is 
lacking. This investigation is concerned with the problem mentioned above \ short discussion is 
devoted to the transitory flow through single capillaries, while the flow minimum involved is qualitatively 
explained. The different behaviour of the flow through micropores is attributed to the influence of the 
pore spectrum, and a new method of determining the “‘mean pore radius” is given. The effect of the 
mean gas pressure and of the te rrp rature on this flow is examined, and the sometimes yreat influence of 
the slip along the walls and the dependence of the mean free path of the gas molecules on temperature 
is elaborated. Finally, the conditions are defined under which surface diffusion may come to play an im 


portant part in the direct flow through micropores. 


Zusammenfassung Die Stromung von Gasen im Ubergangsgebiet zwischen laminarer (PoIseUILLE-) und 
molekularer (KNUDSEN-) StrOmung ist fir die Vakuumtechnik, die in den modernen chemischen Betrieb 
immer starker Eingang findet, von erheblicher Bedeutung Auch treten in diesem Ubergangsgebiet be 

sondere Trenneffekte auf*. Noch liickenhaft sind unsere Kenntnisse tiber die Strémung durch f inporige 
Stoffe, die sich haufig bereits bei Atmospharendruck in diesem Ubergangsgebiet befindet. Die vorliegence 
Untersuchung befaBt sich mit diesem Problem. Die UbergangsstrOmung durch Kinzelkapillaren wird 
kurz besprochen und das dabei auftretende DurchfluBminimum qualitativ verstandlich gemacht. Das anders 
artige Verhalten der StrOmung durch feinporige Stoffe wird auf den Einflu® des Porenspektrums zuriick 

gefuhrt und eine neue Methode zur Bestimmung des ,,mittleren Porenradius’’ angegeben. Die Abhangigkeit 


dieser StroOmung vom mittleren Gasdruck und von der Temperatur wird untersucht und der unter Um- 


* H.-D. Beckey u. W.E. Grorn, Naturwiss. 1951 38 Gasgemische, bei dem sich die leichte Komponente am kalten, 


558. Z. Naturforschg. 1952 7a 474, entdeckten im Uber- die schwere am warmen Kapillarenende anreichert (,,negativer 
gangsgebiet an Kapillaren, die sich in einem Temperatur- Effekt™). Im Bereich der reinen KNupsen-Strémung konnte 
gefalle befanden, einen bisher unbekannten Trenneffekt fir der Effekt zu einem ,,positiven’’ umgekehrt werden 
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Stromung von Gasen durch feimporige Stott 


standen erhebliche EKinfluB der Gleitung sowie der Temperaturabhangigkeit der mittleren freien Weglange 


der Gasmolekeln herausgearbeitet 


Zum SchluB werden die Voraussetzungen angegeben, unter denen eine 


Oberflachendiffusion an der direkten Durchstromung feinporiger Stoffe merklichen Anteil gewinnen kann. 


1. PROBLEMSTELLUNG UND MESSMETHODE 

a) Strémung durch Einzelkapillaren 

Fiir die laminare Strémung von Gasen durch ein 

einzelnes zylindrisches Rohr bzw. eine einzelne Ka- 
pillare gilt bei 1 r nach HaGEen-PorseviLLe 

np ar 


| 
ip 8, L' RT’? 


thy lp Mengenstrom, bezogen auf die Einheit der 
Druckdifferenz Ap: r, L 
Gaszahigkeit: R 


mittlerer Gasdruck: 7 


tadius bzw. Lange der 


allgemeine (Gas- 


Kapillare; 
konstante : p absolute 
Temperatur;: .1 mittlere freie Wegliinge (m. f. W.) 
der Gasmolekeln in der Kapillare beim Druck p. 

Stromung von Gasen 


einen gréBeren Gastransport als nach obiger Beziehung 


Bei der Untersuchung det 
durch Réhren fanden zuerst Kunpr und WarRBURG 


zu erwarten. Die Autoren konnten diese scheinbare 
Verringerung der Gaszahigkeit darauf zuriickfiihren, 
daB an der Wand ein Geschwindigkeitssprung auftritt, 
der durch die endliche GréBe der m.f.W. A verur- 
sacht wird. Unter Beriicksichtigung dieser Gleitung 


laBt sich der Gasdurchsatz genauer erfassen durch die 


Beziehung 
| 4- (2) 
lp Sy L RT 
Hierin bedeutet £ den Gleitungskoeffizienten, der 


etwa von der GréBe der m. f. W. ist: 


const (2a) 


Der numerische Wert der Proportionalitatskonstante 
hingt davon ab, in welcher Weise man die m. f. W. 
aus den experimentellen Ziahigkeitsdaten berechnet 
lval. Gl. (7) und (7a)}. Da die m. f. W. von Gasen 


und Diimpfen bei normaler Raumtemperatur und 


p atm in der GréBenordnung bis em 
liegt, iibersteigt hier der EinfluB der Gleitung einige 
Prozente in Poren vom Radius ; bis 10 
Bei feinporigen Substanzen mul daher unter diesen 
bereits 


auch bei Atmospharen- 


Bedingungen, d.h 


druck. die Gleitung im allgemeinen  beriicksichtigt 


werden 


Ist dagegen die m.f. W. der Gasmolekeln grob 


gegen den Kapillarradius (1 r). so ergibt sich der 
Gasdurchsatz als sog. KxNupsensche Molekularstré- 
(3) 
Ip 3 7 
mit 
2RT 


— mittlerer Molekulargeschwindigkeit 2) 


In diesem Falle gilt fiir die Temperaturabhangigkeit 
des Mengenstroms zweier Gase verschiedener Mol- 
massen und 
| M, 7’, 
Mit H,, O, und CO, hat Knupsen [2] den Gasdurch- 


| T,. (3a) 


satz durch Glaskapillaren in Abhangigkeit von den 
verschiedenen Parametern eingehend untersucht. Bei 
Drucken Giltigkeit der 


Im Gebiet geringerer Drucke stellte er jedoch 


hdheren bestitigte er die 
Gil. (2). 


bei .1 » 3r ein Minimum des Durchsatzes fest. Zu 


Abb. 1. 
Funktion des mittleren Druckes nach KNUDSEN [2]. T ist 


das durch 24 parallele Kapillaren von 2 cm Lange und 


Stromung von CO, durch Glaskapillaren als 


wo C™ Radius unter der Druckdifferenz 1p = 1 dyn em* 
sekundliche flieBende Gasvolumen, letzteres gemessen 
bei p 1 dyn cm? und der Kapillarentemperatur. 
noch kleineren Drucken hin erfolgte wieder ein Anstieg 
um etwa 7°, vgl. Abb. 1. 
Diese MeBergebnisse konnte KNUDSEN mit guter 


Genauigkeit durch eine empirische Formel wieder- 


geben 
= l p 
4 
ap (4) 
mit 
M 
» = é 
c= (4a) 
Hierin sind die Gln. (2) und (3) als Grenzfille ent- 
halten 
Cp : ap h (5) 
ip 
n 
Cp h (Da) 
tp 


Durch Vergleich mit den Gln. (2) und (3) ergibt sich: 
rt l 


3 


(6) 


-|22- (Ga) 
MRI 

In der Form der Gl. (5) kommt besonders deutlich 

zum Ausdruck, daB es sich bei dem Gleitungsglied 

O.S81+ 4 (mit 6 nach G1. (6a)] um einen vom mittleren 


YRS 


Vol 
%, 
~ 
> 
q 
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Gasdruck unabhingigen Zusatztransport handelt, der 
mit steigenden Werten dieses Gasdrucks keineswegs 
verschwindet, sondern lediglich in seinem Verhaltnis 
zum Haupttransport a- jp immer kleiner und schlieb- 
lich vernachlassigbar wird. Setzt man die Gleitungs - 
glieder von Gl. (5) und Gi. (2) einander gleich und be- 


nutzt die Beziehung 
M 


M+, (4) 


A 


l |2 Nao? (7a) 
AY der Partikeln je em?*: a 
) 


StoBdurchmesser), so folgt 


yaskinetischer 


fiir den Zusammenhang 


Schema t psc he Aufteilung einer Stromungs 


in und Knupsen-Ante yur 
Ableitung der GI. (10 


zwischen Gleitungskoeffizient und m.f. W. bei Stré- 


mung durch kreisrunde Réhren [GI (2a) 


(2b) 


Mit dem aus der gaskinetischen Ableitung der GI. (7) 


zu gewinnenden Faktor A '.* hat man also in 
diesem Falle 
1.37-A°%. (2c) 


b) Zum 


stromung von Einzelka pillare ’ 


ly stindnis de Minim hy de Durch. 


Nach der Methode zur Auffindung von Extremwerten 
erhalt man fiir den DurchfluB nach GIL. (4) ein Mini 


mum hei ¢ p Fihrt man in die Definitions. 


Gil. (4a) fiir ¢ die Gaszihigkeit nach GI. (7) ein (mit 
A ',). so folgt fiir die m.f.W.im Durchflub- 
minimum 


310.7 (8) 


Phil 
reine Gase abgeleitete 


Le hrb 


lautet dieser von S. CHAPMAN, lrans 
Roy. Soc. London A 1915 216 279 fiir 
Faktor A 0,499, vel. auch A 
Physik, Bd. Il, 1, 1950, S. 328 
** Fir dic 
gegen 2/3 .A 
loc. cit. S. 387: 


Crenauer 
Chem 
Leipzig 
Platten gilt da- 
nach Gl. (7a)), val \ 
Hand. u 


Leipzig 1939 


Stromung zwischen ebenen 

(mut 
K. F. Herzrecp in 
Physik, Bd. ITI, 2, 


KEN, 
Jahrbuc h der 
Chem 
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Das Minimum in Abb. 1 liegt bei einem mittleren 
0.3 Torr. Hier berechnet sich die 
m. f. W. des CO, nach Gl. (7a) mit ¢ nach Tabelle | 
(OC) zu: Ape 


bei der in Abb. 1 dargestellten Kurve ; 


Gasdruck von j 


9.9 em. Da der Kapillarradius 
3.3° em 
betrug, wird obige Beziehung gut bestitigt. 

Kine strenge Theorie zur Erklirung dieses Durch. 
nicht, Es 


jedoch wenigstens qualitativ verhaltnismaBig leicht 


fluBminimums existiert bisher sich 
verstehen, wenn man beriicksichtigt, daB im Uber- 
yangsyebiet von der KNupsEeNschen zur 
schen Strémung beide Transportarten lings des be- 
treffenden Stroémungsweges hintereinander geschaltet 
Das Druckgefille 


Mengenstrom 


sind lp L fiir einen bestimmten 


nimmt bei 
nach Gl. (1) mit sinkendem Druck j zu. Bei kleinen 
Drucken wird es daher gréBer als das Druckgefiille der 
KNUDSEN-Strémung fiir den gleichen Mengenstrom. 
Wenn 


Drucken ausgehend sich in Teilen des 


das nach Gl. (3) vom Druck unabhangig ist. 
von kleinen 
StrOmungsweges mit steigendem Gasdruck PoIsEvIL- 
Lesche Stroémung anbahnt, sinkt der auf die Kinheit 
des Druckgefilles bezogene Mengenstrom daher 
zunichst unter den fiir reine 
Wert ab. Welcher Anteil 
Stromungsweges Druckabfalls bei 


stimmten mittleren Gasdruck noch reine _Knupsen- 


veltenden des gesamten 


bzw einem be. 
StrOomung enthalt, und welcher bereits zur PoIsEUILLE- 


StrOémung tendiert, hingt wesentlich von dem Ver- 
hiltnis der m.f. W. bei diesem mittleren Gasdruck 


zum Kapillarenradius, ab 


Auf Grund dieser Vorstellung kommt man verhalt- 
nismaBig einfach zu der KNupsenschen Beziehung (4). 
In Abb. La stelle p (7) schematisch den Druckverlauf 
lings des Strémungsweges 0 r L. dar, in dem 
am linken Ende praktisch reine PorseviLLe-Strémung. 
am rechten praktisch reine KNUDsSEN-Strémung herr- 
schen mége. Extrapoliert man den Druckverlauf von 
diesen beiden Enden her zur Mitte hin, so teilt der 
Schnittpunkt A den gesamten Strémungsweg L und 
den gesamten Druckabfall 1p in je zwei Summanden 


L, und Ly baw 


und Ap, vgl. Abb. La. Ver. 
einfachend wird nun angenommen, die PorseviLLE- 
Stroémung erstrecke sich von « — 0 bis zum Punkte A 
(x L,,) heran, von dort ab herrsche dann Knupsen- 
Stromung bis L. Im Teil gilt 
dann fiir den die gesamte Strémungsstrecke 1 durch- 


setzenden Mengenstrom 


a-L-p,: (9) 
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Voll 
No. 6 


mittlerer Gasdruck lings 
, in Abb. la. 
schen Teil gilt fiir denselben Mengenstrom: 


mit a nach Gl. (6) und Pp 


der Strémungsstrecke Im KNUDSEN- 


n=b-L.- (Ma) 
ly 
mit 6 nach Gl. (6a). Da nun L/L L,L 1 ist, 
A 
folgt aus (9) und (9a) 
n= a-p,-A,p+b-Ayp 
oder 
ap: 5. (10) 
Ip Ip Ip 


worin jetzt p den tiber die gesamte Strémungsstrecke L 
gemittelten Gasdruck darstellt 

Wie bereits oben erwihnt, hingt die Einteilung 
des gesamten Druckabfalles in einen PorseEUILLEschen 
KNupseNschen Anteil Verhaltnis 
Ajr ab (1 ist hier die m.f. W. beim Druck j). Da 
1, p/ Ap fiir groBe A/r-Werte verschwindet, fiir kleine 


und einen vom 


sich dem Wert | nahert, wihrend A, p/A p sich gerade 
umgekehrt verhalt, liegt hier der Ansatz 


App l 1 ) Ar 
pt und 


Ip l+aAr Ip 


nahe, in dem « eine empirische Konstante bedeutet. 
Der Faktor p,/p der mittleren Gasdrucke in Gl. (10) 
nimmt mit steigendem A/r (abnehmendem p) von 
Wert | Ahnlich 
Faktor verhilt sich der EinfluB der Gleitung, die 


dem ausgehend zu. wie dieser 


als druckunabhangiger Paralleltransport zur reinen 


PorIsevILLe-Strémung auch bei héheren mittleren 
Gasdrucken grundsatzlich zu beriicksichtigen ist, val. 
Gl. (2). Beide Effekte kénnen gemeinsam durch einen 
Faktor (1 
Gl. (10) beriicksichtigt werden, in dem £ eine 


Konstante darstellt. 


PA/r) bei dem in 
zweite 


Gl. (10) erhalt hiermit die Ge- 


stalt 
- 1+ aAr 
Sie laBt sich leicht auf die Form: 
x a 
l 
= Zz 4 
an+b.- 1] 


bringen, die mit der Knupsenschen Beziehung (4) 
Durch Vergleich hiermit ergeben sich 
0,322 (= 1/3,10) und 
Kurvenverlauf in Abb. 1 
im gesamten Druckbereich gut wiedergegeben werden 
kann. Fiir A r Gl. (10) in die 
Gl. (2 (8 —a)- A. 

Das Minimum des Mengenstroms in Abb. 1 kommt 


identisch ist. 
die Konstantenwerte zu: x 
denen der 


p 5.80, mit 


geht Form der 


liber, mit 


nun dadurch zustande, daB bei geringen Drucken 
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ap<b wird. Dies entspricht dem oben erwahnten 
Fall, daB die Potsevr_Le-Strémung bei kleinen Druk- 
ken fiir einen bestimmten Mengenstrom ein héheres 
Druckgefille 
Dadurch iiberwiegt beziiglich des Druckverlaufs in 
Gl. (11) der Nenner (1 
d.h. 


n/Ap sinkt mit steigendem Druck, 


benétigt als die KnNupsEen-Strémung. 
r/x A), der mit steigendem 


Druck zunimmt: der bezogene Mengenstrom 


von p 0 aus- 
gehend, zunichst unter den Grenzwert 4 ab. bevor 
er mit weiter wachsendem Druck infolge der Zunahme 
des Gliedes a jp wieder ansteigt *. Im DurchfluBmini- 
mum selbst ist a Pp 0,049 5b, wie sich aus den Gln. 
(8) und (11) ergibt. 

Bei der Strémung durch einen Spalt zwischen 
parallelen Platten (Lange groB gegen m. f. W.) tritt 
ebenfalls ein Durchflu@minimum auf, wie W. Gaepe ** 
feststellte, dagegen nicht bei Strémung durch Blenden. 
In letzterem Falle (Strémungsweg kurz gegen m.f. W.) 

. . . . 
kénnen sich im Ubergangsgebiet die beiden rransport- 
arten nicht in der hier angegebenen Weise iiberlagern. 

* Dies fihrt zu der gaskinetischen Vorstellung, daB bei 
den ersten Zusammenst6Ben der Molekeln im 
steigendem Druck, solange deren Zahl noch klein ist gegeniiber 
WandstéBe, die Riickstré 
wiegt, Ahnlich wie bei im Strémungsweg angebrachten Hinder- 
die Zahl der GasstéBe vergleichbar wird 

WandstéBe, die 
zeitig anbahnende gerichtete Bewegung im Sinne einer hydro- 


Gasraum mut 


derjenigen der zunachst iiber- 


nissen. Erst wenn 


mit derjenigen der uberwiegt sich gleich- 
dynamischen Strémung, auf Grund deren dann der Durch- 
flu8 mit wachsendem Druck wieder ansteigt. 

Zusatz hei dir Korrektws 
schienen zwei Arbeiten von J. W 
1952 7a 533, 542, 
Riickstreuung 


Wahrend der Drucklegung er- 
Hisy und M. Pant, Z. Na- 
turforschg in denen der hemmende EinfluB 


unter der vereinfachenden Annahme 


retisch behandelt 


dieser 
Ap Grund einer friiher 
abgeleiteten Streufunktion fir ZweierstéBe, J. W.Hrey und 
M. Pant, Z. Physik 1951 129 517, 130 348 und im AnschluB 
Diffusionstheorie von W.G. Pottarp und R. D. Pre- 
Physic. Rev. 1948 73 762 Fiir kre Kapillaren 
und ebene Spalte ergibt sich die Abnahme des Megenstroms 
Mitteldruck Abb. 1 


in Ubereinstimmung mit den MeBwerten von KNupsEN und 


wird 


an die 


SENT, isrunde 


mit steigendem links des Minimums in 


GAEDE Die Theorie beschriankt sich auf den Druckbs reich, 
in dem die m. f. W. noch groB ist gegen den Kapillarenradius 
(A 20 r) bzw. den Plattenabstand des Spaltes, das eigentliche 


Ubergangsgebiet mit dem Durchsatzminimum wird von ihr 


nicht erfaBt. Fir Lochblenden ergibt sich (unter den Voraus- 


setzungen Ap<p, A r) in Ubereinstimmung mit der Er- 
fahrung eine stetige Zunahme des Durchsatzes mit steigendem 
Mitteldruck, was auch noch fir kurze Kapillarenstiicke gilt; 
oberhalb L = 8r liefert die 


verminderung mit steigendem jf. Im 


Theorie Durchsatz- 


AnschluB hieran fiihren 


erst eine 
die Autoren das Fehlen eines Durchsatzminimums bei der Durch- 
stromung poréser Medien darauf zuriick, daB in diesem Falle 
AL Sr gelten dirfte, wenn man unter AL die geradlinige 
Strecke eines Porenstiickes zwischen Knicken versteht. 

** W. Gaepe, Ann. 1913 (4) 41 289. Der Autor 


vermutete als Ursache des DurchfluBminimums eine an den 


zZwel 


Phys. 


Platten haftende Gashaut. 


@ 
* 
OL. 
95 2? 
= 


and W. 


Wick? 


c) Grundsitzliches zur Strémung durch porose Sub- 


stanzen 

Wesentlich schwieriger ist der Gastransport durch 
porése Substanzen zu erfassen. Es liegt in der Struk 
tur dieser Kérper, daB die Porenradien und -langen 


oft um Zehnerpotenzen streuen Immerhin abt sich 


aus Messungen des Gasdurchflusses bei verhaltnis 
mabBig hohen und bei sehr kleinen mittleren Drucken 
ein fiir Zwecke Mittelwert des 


Porenradius erhalten. Mibt man bei zwei Werten p, 


manche niitzlicher 


und p, héherer Gasdrucke den Durchsatz #, und n”, 
bei gleicher Druckdifferenz 1 p, was einer Bestimmung 
der Neigung der nach rechts auslaufenden Geraden 
in Abb. 1 entspricht, so la4Bt sich durch zweimalige 
Anwendung der Gl. (2) der Gleitungsanteil in erster 


Naherung eliminieren: 


(Z, Zah!l der Poren mit dem Radius r; und der 
Linge L,.) 

Kombiniert man diese beiden Messungen mit einer 
dritten, die bei sehr kleinen Drucken im Giltigkeits- 
Gl. (3) (Druckdiffe- 


1, p, Mengenstrom so erhilt man den er- 


bereich der durchzufiihren ist 
renz 


wihnten Mittelwert des Porenradius zu 


A 
> & 
4i° 7 
| 
(12) 
16 “ Ny p 
3 Py Ps Ip "| 


Bei mittleren Drucken tritt im allgemeinen lamin- 
are SchichtenstrOmung nach Glei- 
tung, Knupsensche Molekularstr6mung und gegebe- 
den Porenwan- 


nenfalls Oberflachendiffusion an 


dungen zum Teil nebeneinander auf, zum Teil gehen 
diese verschiedenen Transportarten in verschiedenen 
Was hierbei 


als héherer, mittlerer oder kleiner Gasdruck anszu- 


Abschnitten einer Pore ineinander iiber. 
sprechen ist, haingt weitgehend von dem speziellen 


Charakter des betreffenden Porengefiiges ab. Es 
kommt hier vor allem auf die GréBenordnung des Ver- 
haltnisses der m. f. W. zum mittleren Porenradius an, 
Gl. (12) 


anderen geeigneten Mittelwertsbildung erhalten kann. 


wobei man letzteren etwa aus oder einer 


In der vorliegenden Arbeit wurden daher die Stré- 
mungsverhaltnisse in porésen Substanzen nach folgen- 
den Punkten untersucht: 


1. Anteil und EinfluB der Gleitung am Gesamt- 


transport bei A < 7. 
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2. Druckabhangigkeit des Gesamttransportes bei 
kleinen Drucken (1 > 7). 

3. Einflub der Temperaturabhangigkeit der m.f.W. 
Anteil 


Stromungs- 


$. Voraussetzungen fiir einen merkbaren 


der Oberflichendiffusion am gesamten 


transport 


d) A pparatur und Versuchsfiihrung 

Die zur Klaérung dieser Fragen benutzte Apparatur 
war verhaltnismaBig einfach, vgl. Abb. 2. Sie bestand 
im wesentlichen aus zwei auswechselbaren Glaskolben 
A, und A,, die tiber ein Leitungssystem und das 


StrémungsgefaB S miteinander verbunden waren. 


Verschiedene Vorratsbehalter V enthielten die vor- 
gereinigten Versuchsgase. Als porése Substanz diente 
einmal eine Jenaer Glasfritte der Sorte G4 (4,1 mm 

Hilfe 


von GI. (12) zu 2.2- 10-* em bestimmt wurde (Katalog- 


10> em). 


2mm stark), deren mittlerer Porenradius mit 


angaben 4,5 bis 7.5- Bei den Messungen 
mit diesem Filter betrugen die Volumina der Kolben 
A, und A, je etwa 61. Eine zweite porése Probe be- 
stand aus feinteiligem Nickelpulver, das aus Nickel- 
oxyd durch vorsichtige Reduktion mit H, bei 280°C 
gewonnen und mit 5000 atm zu einem kleinen Zylinder 
von 8 mm und 8,5 mm Héhe gepreBt worden war. 
Dieser PreBzylinder mit einer scheinbaren Dichte von 
5.0 g/em® wurde in das Ende eines Glasrohres ein- 
gekittet (an der Stelle F in Abb. 2); fiir den mittleren 
geschilderten 


3.0- 10-4 em. 


Porenradius ergab sich nach der oben 
Methode mit Hilfe Gl. (12): 7 
In diesem Falle war es zweckmiabig, die Kolben A, 


von 


und A, der MeBapparatur (Abb. 2) nur je einen Liter 
groB zu wihlen. 

Die Geschwindigkeit der GasstrOmung durch die 
porésen Fritten bzw. Filter wurde durch Verfolgung 
des Druckausgleichs zwischen den Kolben A, und A, 
Eine Mano- 
eingestellte Druckdifferenz 1 p, 


instationér gemessen. anfinglich am 
meter D (Abb. 2) 
setzt eine Ausgleichsstrémung in Gang, die mit ab- 
nehmender Druckdifferenz dieser proportional ab- 
klingt : 

A l 


l d(Ap) 
RT 


n Ip f(p, T) 
Hierin bedeuten V, und V, die Volumina der Kolben A, 
tohrleitungen bis zur 
Die GréBe T) 


der be- 


und A, einschlieBlich der 
Strémungsfritte bei F in Abb. 2. 
Charakter 


treffenden 


eines Strémungsleit wertes 


Fritte ist 


vom 
porésen wihrend einer Einzel- 
messung des Ausgleichsvorganges natiirlich konstant, 


im tibrigen abhangig von der Art der Strémungsfritte 


2Sb 


q 
a 
> 
“ps 
= 
Ip 84 RT 
i 
j 
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und des Gases sowie vom mittleren Gasdruck und 
der Temperatur innerhalb des durchstrémten Poren- 
Druck- 
In Ap(t)/Ip, gegen 


gefiiges. Wird die zur Zeit ¢ gemessene 
differenz Ap(t) in der Form 
die Zeit ¢ aufgetragen, so erhilt man eine Gerade. 
deren Neigung den Wert von 1/4 p- d(A p)/dt angibt, 
der zur Ermittlung des auf Ip | Torr bezogenen 
Mengenstroms #/Ap in Gl. (13) einzusetzen ist. 


Abb. 2. Stromungsapparatur. V Gasvorratsbehalter; 
A,, A MeBkolben; M Manometer fiir Absolutdrucke: 
S Stromungsgefa8 mit F poréser Probe (Glasfilterfritte 
oder eingekitteter PreBzylinder aus Ni-Pulver); K Kuhl 
fallen; D Differentialmanometer. 


2. MESSERGEBNISSE UND AUSWERTUNG 
a) EinfluB der Gleitung auf den Gesamttransport bei 
A<? 
Die mit dem Ni-PreBzylinder und dem Jenaer Fil- 
ter G4 bei 0°C mit H, bzw. N, zwischen Atmo- 
spharendruck und 2,5 Torr gemessenen Gasdurchsitze 
sind in Abb. 3 und 4 dargestellt, und zwar die GréBe 
nj Ap = f(p, T), Gl. (13), gegen den mittleren Gas- 
druck p. Nach der Hacen-Potsevitieschen GI. (1) 
sollte die Extrapolation der bei héheren 
Drucken gemessenen Gerade auf p—> 0 durch 


den Koordinatenanfangspunkt gehen. Das 


samttransportes bei Atmospharendruck erreichen 
kann, 

Besonders interessant ist in diesem Zusammenhang 
eine Untersuchung von W. Brérz und H. Srencier 


3}, in der die Strémung einer Reihe von Gasen durch 


eine Probe eines feinporigen technischen Fiscner- 
Tropscu-Katalysators bei 350 Torr gemessen 
wurde. Tabelle 1 zeigt in den Spalten 1—4 die unter 


‘ 


Stromung von Wasserstoff bei 0° C durch den Ni- 
PreBzylinder (7 3,0 em). 


Abb. 3. 


Abb. 4. Strémung von Stickstoff bei 0° C durch die 


Glasfritte (¥ 2,2-10-*em). 


Anwendung von Gl. (1) fiir reine Portsevmwe-Stré- 


mung gewonnenen Ergebnisse. Die zuniachst uner- 


wartete Streuung der Werte des Ausdrucks y 


ip. 


Tabelle 1. EinfluB der Gleitung auf die Strémung verschiedener Gase durch 


eine por OR Katal ysator probe hes p 350 Torr nach Me ssungen von 


BroTz und SPENGLER 


AusmaB der Abweichung von dieser Be- 


dingung zeigt den erheblichen EinfluB der : 
Gleitung an, die nach Gl. (2) fiir eine Beschrei- 
bung der experimentellen Werte herangezogen Megas 
werden muB. 

Das Wort Gleitungs-,,Korrektur’, das 
man haufig in der Literatur findet, sei hier ab- CH, 
sichtlich vermieden. Wie bereits oben betont, N, 
handelt es sich bei Porenradien von 10-4 em NH, 
und kleiner keineswegs um eine Korrektur- CO, 
gréBe, sondern um einen Strémungsanteil, H, 

CO 


der durchaus die GréBenordnung des Ge- 


2 3 f 5 6 
mol Poise Ip r = 
l 4 


sec dy nycm-* 


2,34 108 2,52 2,38 1,06 
1,55 175 2,7 2.69 1.01 
2,56 100 2,56 2,24 1,14 
1,52 145 2,20 2,11 1,04 
4,82 SS 4,23 4,18 1,01 
1,67 177 2,96 2,7 1,09 


a | 
A 
| 
Ihe 
| 
LJ 
J 
OL. 
| 
| | | 
952 «109 4 
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in Spalte 4, der nach Gl. (1) fiir jedes MeBgas derselbe 


sein sollte, erweist sich als spezifischer EinfluB der 


Beriicksichtigt 


Spalte 5 angegebenen Gleitungsfaktor nach Gl. (2), 


Gleitung. man diese durch den in 


so ergeben sich die in Spalte 6 verzeichneten Werte, 
die nun in der Tat als von der Natur des MeBgases un- 
abhangig angesehen werden kénnen 

Zur Berechnung des Gleitungsfaktors wurde hier- 
bei ein mittlerer Porenradius von 7 450+ 10-5 em 


gewaihit, den man aus den Angaben der Autoren 


Fiir den Glei- 
nach Gl. (>< 1-37 


gesetzt. Die Werte fiir die m.f. W. bei 20° C und jp 


durch Anwendung von GI. (12) erhailt. 
tungskoeffizienten wurde 
350 Torr wurden aus den in Spalte 3 angegebenen 


Zihigkeiten nach Gl. (7) berechnet. 


b) Druckabhangigke if de samttrans ports hei kl ine a 
Drucken (A > 7) 


Der experimentelle Befund fiir die untersuchten Pro- 


ben ist bereits in den Abb. 3 und 4 enthalten. Die 
010% 
A=7 
20 600 
Abb. 5. Strémung von Waaserstoff bei 20° C durch einen 


PreBling aus Al,O,-Pulvet 2-10° en ach 


CARMAN [4 


dort eingezeichneten Pfeile geben an, bei welchem 
Druck die m. f. W. gleich ist dem nach GI. (12) be- 
Bei dem Ni-Preb- 


zvlinder nimmt in der Umgebung dieses Punktes die 


rechneten mittleren Porenradius 
Neigung der gemessenen Kurve deutlich ab, ein Durch- 
flubminimum wie bei einzelnen durchstrémten Kapil- 
laren (Abb. 1) wurde jedoch in keinem Falle beob- 
achtet. P.C 
liber die Strémung von H, durch PreBlinge aus fein- 


teiligem Al,O,-Pulver fest, daB die im Gebiet héherer 


CaRMAN [4] stellte bei Untersuchungen 


Drucke gemessene Gerade in einzelnen Fallen die 
MeBpunkte bis zu niedrigen Drucken (1 r) herunter 
ohne merkbares Abknicken wiedergibt, val. Abb. 5. 

Bei porésen Stoffen bietet somit die Proportionali- 
tat des Strémungsdurchsatzes zum mittleren Druck 


keine Gewihr dafiir, daB es sich um reine Porsevr..e- 


Strémung handelt. Vielmehr kann hier durchaus eine 
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scheinbare Druckabhaingigkeit der KNupsEN-Stré- 
mung ohne Minimum vorliegen, hervorgerufen durch 
die Uberlagerung der verschiedenen hintereinander- 
und parallelgeschalteten Porenweiten innerhalb der 
betreffenden Strémungsfritte. Uberlagert man z. B. 
mit Hilfe der KNupsenschen Beziehung (4) die Stré- 
mung durch zwei zylindrische, parallel verlaufende 


Zehner- 


potenz unterscheiden, so kann im Gesamtdurchsatz 


Porensorten, die sich im Radius um eine 


das bei der engeren Porensorte auftretende Minimum 


verschwinden,. vgl. Abb. 6. 


Einflup der m pr raturabhdngigk it der m.f. W. 


Glasfritte Gd Die wurden bei einem 


mittleren Druck von 3 Torr durchgefiihrt, dem eine 


Messungen 


/ 4 


Abb 6 


denen 


Bei der Durchstromung pordser Substanzen, in 
Poren verschiedener Weite einander parallel ge- 


schaltet sind, kann das von KNUDSEN in einzelnen Ka 


pillaren bei 3 


schwindet 


festgestellte DurchfluBminimum ver- 


Kurve a: Durehflu8 je Einheit des Druck 


gefall von H, bei cinem Gesamtdruck von 760 Torr 
und © C durch 1000 Kapillaren vom Radius 107-4 «em 
und der Lange J lem. Kurve } DurchfluB durch 


denselben Be 
Parallel 
mittels Gl. (4) 


eine Kapillare vom Radius 107% em unter 


dingungen. Kurve a 6: GesamtdurchfluB bei 


schaltung beider Porenarten (berechnet 


m.f. W. von | bis 2- 
entsprach. Tabelle 2 enthalt die Ergebnisse der Ver- 
bei 


em der verwendeten Gase 


suche 


Tahelle Clasatrom ung dure Glastritt C und p 3Torr 


(n, Ip in cm® NTP) sec Torr) 


ny Ip: ny Ip-| M 
He 7.4 
N. 04 OLOLOSO 
CH, 6S O.OL071 
C Hy 1,46 
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Die Werte der ersten drei Gase passen sich im zu- 
liissigen Fehlerbereich der | M-Beziehung nach GI. (3a) 
an. Aus dem Mittelwert von (n/4 p)- | M = 0,01073 
ergibt sich dagegen fiir Butan ein UberschuBtransport 
von 3,8°%. Wicker und Vorer* fanden bei diesem Gas 
unter Verwendung einer Glasfritte G5 einen Uber- 
schuBtransport von 5,3°% bei 0° C. 

Die Messungen wurden sodann auf ein gréberes 
Temperaturintervall ausgedehnt, um die Abhingig- 
keit von dem Faktor | MT in Gl. (3a) naher zu pri- 
fen. Wie die in Abb.7 eingetragenen MeBpunkte zeigen, 
wurde auch bei N, und CH, ein mit sinkender Tempe- 
ratur ansteigender Chers« huBtransport im Sinne einer 
Abweichung von der | MT-Beziehung festgestellt. 


berechnet gemes: 


780 


Abb. 7. 


durch 


Temperaturabhangigkeit des Gastransportes 
bezogen auf die | M7’-Beziehung 


3 Torr. 


die Glasfritte, 


(3a) bei p 


Ni-PreBzylinder — Bei dem gepreBten Ni-Pulver 
betrug die m. f. W. der MeBgase bei 0° C 
wihlten mittleren Gasdruck von 3 Torr ebenfalls noch 
etwa Zehnfache 
radius (3-10°4em). In 


MeBpunkte eingetragen. 


und dem ge- 


durchschnittlichen Poren- 


Abb. 8 sind 


das des 
die erhaltenen 


Hiernach ist bei tiefen 


Temperaturen das Produkt n/4p. | MT auch in 


diesem Falle nicht konstant, sondern steigt zum Teil 
erheblich an, beim CH, z. B. bei der Temperatur 
fliissiger Luft um etwa 50% gegeniiber 0° C, beim H, 
um etwa 

Diese Zunahme des Strémungsdurchsatzes mit 
sinkender Temperatur laBt sich in allen hier gemessenen 
Fillen quantitativ wiedergeben, wenn man die Tem- 
peraturabhiangigkeit der m. f. W. der betreffenden Ga- 
se beriicksichtigt. Man stellt sich hierzu zweckmaBig 
die porése Fritte in vereinfachender Weise als ein 

* E. Wicke und H. Voter, Angew. Chem. (B) 1947 19 94. 
Dort wurde dieser Zusatztransport fir Oberflachendiffusion 


gehalten. Wie die nachfolgende Diskussion zeigt, beruht er 
dagegen tatsichlich auf einer Verschiebung zwischen Knup- 
sEN- und Porsevrtie-Anteil, hervorgerufen durch die An- 


derung der m. f. W. mit der Temperatur. 


22 Chem. Eng. Sei. Vol. 1 


Stromung von Gasen durch feinporige Stoffe 


System von Z parallel geschalteten Poren vom gleichen 

tadius r und der Lange Z vor und wendet auf dieses 
System die vollstindige KnupsEeNnsche Beziehung (4) 
bzw. Gl. (1) an. Zur unmittelbaren Wiedergabe der 
MeBwerte bringt man die Gleichung zweckmiaBiger- 
weise auf die Form: 


| MT-n Ip =520-Z2P L 


(14) 


O0,81-kr/A \ 
l+kr/A 


in der die Zahlenwerte so gewahlt sind, daB n/A p den 


Gasdurchsatz in cm’ NTP je Torr Druckdifferenz und 


Abb. 8. 


durch den 


Temperaturabhangigkeit des 
Ni-PreBzylinder, 
Beziehung (3a) bei p = 3 Torr. 


Gastransportes 


bezogen auf die 


Sekunde bedeutet. Die Gleichung enthalt zwei von 
der Gasart und der Temperatur unabhingige Kon- 


Faktor 


r, die den EinfluB des ausschlaggebenden 


stanten: den vor der Klammer und 
die GréBe k- 
Ausdrucks in der Klammer wesentlich bestimmt. 
Durch Anpassung an die in Abb. 7 und 8 angegebenen 
MeBpunkte ergaben sich fiir diese Konstanten bei der 
Glasfritte G4 bzw. dem Ni-PreBzylinder die empiri- 
Werte: ZrL 17+ em? 1,73 

10° em?; ker 10-10°% bzw. 2.3- 10-%em. Die 
m. f. W. StoB- 
querschnitte der Gasmolekeln wurden dem LAaNDOLt- 


Bd. 1/1 


Tabelle 3 sind die gaskinetischen 


schen bzw. 


zur Berechnung der erforderlichen 


Bornstein, 6. Aufl., (1950) entnommen. In 
Durchmesser der 
betreffenden Gase zusammengestellt, ferner die Svu- 
THERLAND-Konstanten fiir den interessierenden Tem- 
und 273 


peraturbereich zwischen 90 abs. sowie die 


daraus berechneten m. f. W. 

Die ausgezogenen Kurven in den Abb. 7 und 8 sind 
mit den oben angegebenen Werten nach Gl. (14) be- 
worden. 


rechnet Der mit sinkender Temperatur 


. 
ansteigende ,,UberschuBtransport“ der Gase wird 
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‘ 
k °7 A 21 
‘ 
| 
| berechnet gemesse “4 
N 
> \ 
9 5 ? > > ~ 
: 


E. Wicker and W. 


Gaskinetische Durchmesser (in 10°° cm). SutTuer 


Tabelle 3 
LAND-Konstanten und mittlere freie Weglangen bei ; 3 Torr 


fin tom) 


273) (572) C°K A(90) 


3.06 36,4 
4.65 91.8 
5,31 118.2 
4.70 27,2 


6,92! 280,08 


T 27: 
273 und 393° K 


hierdurch sehr gut wiedergegeben und eindeutig ge- 


klart 


so groB, daB die Beziehung 


Bei héheren Temperaturen sind die m.f. W. noch 
MT 
Niaherung gilt (fiir Butan erst oberhalb des hier itiber- 


const mit 


= — 


200 


Abb u 


hangigkeit vom 


Stromung von CF,Cl, durch Linde-Silica in Ab- 
Druck bei 
MALHeERBE loc. cit 


Druckdifferenz in Torr; 


nach CARMAN und 


Gasstrom in millimol/sec; 
L Probendicke in em; A Proben- 


querse hnitt in em?. 


schrittenen Temperaturbereichs zu erwarten). Infolge 


der erheblichen Zunahme der Wirkungsquerschnitte 


und Abnahme der m. f. W. nimmt der Gastransport 


mit sinkender Temperatur jedoch gerade in den Poren 
mit gréBeren Querschnitten immer mehr den Cha- 
rakter einer PorseurLie-Strémung an, die den in 
Abb. 7 und 8 gezeigten Anstieg des gesamten Gas- 
durchsatzes hervorruft. 


d) Voraussetzungen fii 


Anteil der 


Oberflichendiffusion am gesamten Striim ungstransport 


In allen vorangehend besprochenen MeBergebnissen 
konnten die Abweichungen von dem erwarteten Stré- 
durch Uberlegungen 
Ein EinfluB von Ober- 
flachendiffusion trat bei der direkten, von einer Druck- 


mungstransport gaskinetische 


befriedigend erfaBt werden. 


differenz hervorgerufenen Gasstrémung nicht in meb- 


Stromung von Gasen durch feinporige Stoffe 


guter 
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barer GréBenordnung auf. Im Zusammenhang hier 
MAL- 
HERBE [5] von Interesse, in denen Oberflichendiffusion 


unter 


mit sind Untersuchungen von Carman und 


Anwendung der Strémungsmethode nachge- 
wiesen werden konnte. Als bestes Beispiel seien die 
Linde-Silica IT bei 


33° C angefiihrt, vgl. Abb. 9 (innere Oberfliche etwa 


Messungen dieser Autoren an 
600m? cm*, mittlerer Porenradius etwa em). 
Als Gas wurde CF,Cl, verwendet, dessen Sitti- 
gungsdruck bei der Versuchstemperatur 660 Torr be- 
Die Kurve a in Abb. 9 gibt den gemessenen 
Durchsatz des CF,Cl, an. Die Kurve b, die den Anteil 


der reinen KNupseN-Strémung darstellt. beruht auf 


tragt. 


Stromungsmessungen mit Luft, die mit Hilfe der Mol- 
massen nach Gl. (3a) auf CF,Cl, umgerechnet wurden. 
Die Differenz zwischen a und 4 interpretierte Car- 
MAN als Oberflichendiffusion, was qualitativ sicher 
zutreffen dirfte. In quantitativer Hinsicht ist aller- 
dings zu bemerken, daB ein Teil des mit steigendem 
Druck auftretenden Zusatztransportes auf ahnliche 
kénnte wie der in 
Abb. 7 und 8 gezeigte Transportanstieg mit sinkender 
Temperatur. Mit steigendem Druck sinkt die m. f. W. 


Ursachen zuriickzufiihren sein 


ab und damit nimmt ein Teil des Gastransportes, ins- 
besondere in den weiteren Poren, den Charakter von 
PoIsevILLE-Strémung an, was i.a. 
Es liegt 
zunehmen, daB der steile Ansteig in Abb. 9 zum Teil 
Einfliisse Dies 


wiirde zugleich erklaren, weshalb der gesamte Zusatz- 


bereits einen Zu- 


satztransport hervorruft. daher nahe, an- 


durch dieser Art verursacht wird. 
transport schneller als proportional dem Oberflachen- 
Konzentrationsgradienten ansteigt, was CARMAN als 
unerwartetes Ergebnis besonders betont. 

Wie aus diesen Versuchen folgt, kann die Ober- 
flachendiffusion bei der direkten Durchstrémung eines 
pordsen Kérpers dann einen meBbaren Transport - 
wenn die innere Oberfliche der 


anteil liefern, von 


GréBenordnung 100 poroser Substanz oder 
grOber ist, was mittleren Porenradien von der GréBen- 
ordnung 10°*em oder darunter entspricht. Ferner 
ist es hierfiir giinstig, wenn der mittlere Druck in 
der Nahe des Siattigungsdruckes des Versuchsgases 
liegt, damit die innere Oberfliche eine hohe Adsorp- 
tionsbelegungsdichte, méglichst in mehreren Molekel- 
lagen, aufweist. Dann wird der dem Druckgradienten 
im Gasraum entsprechende Gradient der Oberfliachen- 
konzentration sehr hoch, und es kann unter Um- 
stinden eine fliissigkeitsihnliche Adsorptionsschicht 
gleichsam ,,flieBen“, was den Gesamtdurchsatz dann 


wesentlich erhdht. 
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0(90) o( A (273) (572) 
H, 7.45 28.8 66.5 
N, 3,23 14,8 36,0 1. 
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Ar 3.15 15.85 30.7 
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Solche Verhaltnisse lagen bei den hier untersuch- 
ten Glasfritten und Ni-PreBlingen in keinem Falle vor. 


Der eine von uns dankt der Ruhrchemie AG., Oberhausen- 
Holten, fiir eine finanzielle Beihilfe. die ihm die Durchfiithrung 


voranstehender Untersuchungen im Rahmen seiner Disser- 


tation ermogli« hte. 


ZUSAMMENSTELLUNG DER FORMELZEICHEN 

= SUTHERLAND-Konstante ° K 

i = Index fiir verschiedene Porensorten 

Index fiir KNupsEen-Strémung 

L == Lange der Strémungsstrecke: em 

M Molmasse ; 

Gasmengenstrom; mol ‘sec 

(lem* NTP see 1/22415 mol sec) 

N Partikelzahldichte; 
Druck, mittlerer Druck, Druckdifferenz: 
1333 dyn/em?) 


p = Index fiir Potsevitie-Strémung 


g mol 


dyn/em*® (1 Torr 


r, 7 Kapillarenradius, mittlerer Porenradius: em 


Book 


Watter G. Bert (Editor); Physical Methods in Chemical 
Analysis, Vol. I. Academic Press Ine.. Publishers, New York 
1950. 664 pp. Price $ 12-00 

La public ation d'une en ve lopédie sur les méthodes phy- 


sIques applicables a analy se chimique repond a un besoin 


certain, mais la rédaction d’une pareille ency: lopédie risque 


d’échouer devant l’énormité de la ta he, en se rdant dans 


les détails, ou de ne donner que des notions pratiquement 


inutilisables sur les méthodes phy siques s’appliquant l’ana- 


ly se chimique moderne. Les douze chapitre s qui torment le 


volume I répondent d'une facon trés satisfaisante lobjet 


d'une pareille en yelopédie et fournissent une revue a jour 


des progrés actuels sans cher: her A rivaliser avec des mono- 
graphies consacrées A chacun des sujets. 


CLARK, W.L 
Davipson, J. HowsMon traitent des applications des rayons X 


Les trois premiers chapitres dis a G. L 


soit aux phénoménes d’absorption, soit aux phénomeénes de 
Nous 


avons particuliérement apprécié les pages traitant de lappli- 


diffrac tion, dans le cas des poudres, metaux ou fibres 


cation des rayons X a l'étude des propriétés physiques (duretés 
proprietes magnetiques et déformation). 

Un chapitre sur la diffraction électronique par L. O. 
trés clair. 


BROCKWAY est plus bref, tout en restant 


Les chapitres sur l’analyse spectrale avee ses diverses 


branches sper trophotographie et colorimétrie par W.R. 
Brope, spectrographie d’émission par J. SHERMAN, spectro 
scopie infra-rouge par H. H. Nietsen et Ropert A. N, 


spectres Raman par J. H. condensent d'une facgon 


remarquablement heureuse l'ensemble de ce qu'il faut savoir 


our AD uer le spectrogra rhe a l’immense Variete des 
pq I I 


problémes analytiques qu’ il permet de résoudre. 
Il reste évidemment moins de place pour des chapitres 
plus classiques tels que les applications de la lumiére polarisée 
par C.D. West et de la réfra tometrie par L. W. Tu TON 
CODSACTCES a Ces 


J.K.Taytor. Mais les 110 pages qui sont 


deux branches de optique, appliquées depuis longtemps 


& l'analyse chimique, sont excellentes. 
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8.315 +107 


allgemeine Gaskonstante 
dyn-+em 
mol - K 
t Zeit; sec 
T = absolute Temperatur; ° K 
mittlere Molekulargeschwindigkeit: cm/sec 


Z Porenzah| 


a, p Konstanten 
Caszahigkeit : Poise 
cm 8e¢ 


Gleitungskoeffizient (Gl. 2): em 


A mittlere freie Weglinge [Gl. (7a)]; em 
o — gaskinetischer StoBdurchmesser (Gl. (7a)]; 


em 
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reviews 
R. D. HEmENREICH consacre 50 pages & la microscopie 
électronique. Elles permettront de s erendre compte des 
ressources de ce nouveau moven d’investigation. Mais il 


serait déraisonnable Wespérer qu’elles suffisent & former les 


operateurs auxquels l’on ¢ onfierait un microscope elec tronique. 
Nous ferons la méme remarque sur les chapitres consacrés 


H.W. Wasusurn. Ce 


chapitre donne une idée parfaitement précise de 


a la spectrometrie de masse par 
ce que 
peut attendre de la spectrometrie de masse pour l’analy se des 
substances « omplexes, 


Le premmer volume de la collection doit étre particuliére- 


ment apprecie par les directeurs de recherches qui tiennent 


a etre tenus au courant des ressources que peuvent leur 


apporter les méthodes les plus modernes de la physique pour 


la solution des problémes analytiques. J. CATHALA 


W. Josr; Diffusion in Solids, Liquids, Gases. Academic Press 

Inc., Publishers, New York, 1952. 558 pp. Price § 12-00. 
Many readers will remember an earlier book by the same 

author under the title “ Diffusion und ¢ ‘hemische Reaktionen in 


festen Stoffen”, which appeared in 1937. Rather than prepare 


a second edition of this, Professor Josr has preferred to 
treat diffusion processes on a broader basis and to neglect the 
subject of chemical reaction in solids entirely, except for 
such cases that could be described in terms of diffusion pro- 
cesses”’. The result is the present volume. 
The first chapter is on the fundamental laws of diffusion 
useful collection of solu- 
the 


This includes 


and contains an extensive and wv ry 
the Fix k 


experimental conditions met with in practice 


tions to diffusion equation covering most of 


solutions not previously published for diffusion in systems 


containing more than one phase, (in the me tallographic sense ) 


a case of considerable importance in industrial practice. 


The se ond chapter describes and gives the theory of the 


various types of structural defects in solids, particularly 


those which are responsible for diffusion and onic conductivity 
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in solids. This is followed by one on the general theory of 
diffusion and conductivity processes themselves 


rhe next two chapters deal at some length with electro- 


lytic conduction and diffusion in ionic crystals and with 
ditfusion in metals and non-polar crystals respectively 
Accounts are given of the experimental methods of investi- 


gation and, as in all the chapters of the book dealing with the 
more experimental aspects of the subject, a very large number 
of numerical results are quoted 

A chapter on solubility, particularly of gases in solids, 
paves the way for Chapter VII, which is on the permeation 
and diffusion of gases in solids 

Chapter VIII deals with the electrolytic mobility of ions 
in solid and molten metals and alloys 

Chapter LX is devoted to surface reactions of metals which 
involve diffusion such as the growth of oxide and sulphide 
films and chemical reactions between powdered materials 
A number of such reactions are selected for detailed treat 


ment, the experimental methods which have be« mploved 


in their study being described and the results obtained dis 


issed in terms of WaGNeER’'s theory, of which a good account 


The next two hapters ure respectively on 
gases and diffusion in liquids. Both give excellent surveys 
of the exp rimental methods and useful collections tf experi 
mental results. The theoretical treatment of gaseous diffusion 
8 restricted to the quoting of a number of formu isetul 
n practical application, for excellent monographs are already 
available on the « ibject, but a fuller account is given of the 
generally less successful attempts that have been made at 


a theoretical understanding of | quid diffusion 
The book concludes with a useful and interest hapter 


on the theory and « xperumental methods of thermal diffusion 


in gaseous and condensed phase mixtures 

Extensive bibliographies and references are given at th 
mls of wl ipter 

Phe emphasis throughout is on the experimental methods 
of diffusion stuches and their results and the book is notable for 
tl irg mber of numencal results which ar juoted: 
for this reason it forms a most useful source of diffusion data 
But in many es the presentation of data is u ritical, 


g made to assess their reliability or a iracy, 


even when disparent results by different investigators of the 


same system are quoted side by sick 

The summarn that are given of experunental methods 
re very satisfactor and mprehensive, and the same is 
enerally true of the theoretical sections of the book There 
ire, however, a number of omissions which are to | gretted 

a book of this size Thus, for example, practica thing 
* sani of surtace liffumon, or of grain bouncdas 
me ti important contribution it may make to a baerved 
diffuson flow in polyerystalline materials at low temperatures 
Dif fuss m cleetrolyt solutions mil receives 
ver scanty attention, and it is unfortunate that so litth 
SL = oh ted to the Kirkendall effect and ew leas to 


DARKEN's analysis of it for these have been central themes 
in recent work on diffusion, especially in metals 

However, in spite of these shortcomings the author has 
largely succeeded in his aim “to assist in planning, evaluating 
and understanding diffusion « xperiments, at the same time 


giving a survey of results to date”. Although ‘to date” is, 
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unavoidably, nearer to 1949 than the present time, references 
are given in an appendix to papers published as late as 195] 

The printing and binding call for no special comment 
except for a misprint” on page 306, the whole of which 
seems to have been printed from ty pe set up for quite another 
book! The extensive reader may find the missing page in 
other publications of the Academic Press, but others will 
find its omission detracts noti« eably from the value of Chap 
ter VII. 

The book may be safely recommended to all interested 
and working in this field as a very useful addition to their 
shelves, but many may feel that its merits, many though 


they be, are hardly adequate to justify its very high price 


A.D. Le CLarre 


A. R. Usseroupe; An Introduction to Modern Thermo- 
dynamical Principles. 2nd Edition, Oxford University Press 
1952 pages 

In this book Professor Unpne LOHDE has given a short but 


statistical thermodynamics, It will be of real value to chemists. 


comprehensive survey of the main fields of application « 


engineers and physicists who wish to know the kinds of problem 
which may be solved by the use of statistical thermodynamics. 
as a preliminary to studying the more complete treatises 
\ great deal of ground is covered at a fair pace and almost 
inevitably there are a few errors and obscurities. For teaching 
purposes the book could be warmly recommended to those 
students who have already taken a course of systematic 
lectures on the foundations of the subject 

The first three chapters are concerned with thermodynami 
functions and the fourth and fifth with the Nernst heat theorem 
In Chapter VI there is an elementary account of the statistical 
mechanics of a system of independent particles, using the 
method of the most probable distribution and Sriruma’s 
theorem. Chapter VII presents a useful survey of the statistical 
thermodynamics of the solid state, including the Desye fre 
quency spectrum, co-operative effects and the Brace and 
WitiiaMs theory of order-disorder transitions in alloys Chap 
ter VIII is concerned with the methods of attaining very low 
tripe ratures 

In Chapter IX there is a short outline of Einstein -Bost 
and Frerei-Dirac statistics. This chapter does not connect 
up very well with the earlier discussion in Chapter VI of the 
BOLTZMANN statistics and in particular there is no discussion 


of the important factor log n! in the formulae for fre energy 
and entropy However it is notoriously difficult to carry 
out a really satisfactory treatment, at an elementary level, 
of the quantum statistics of non-localised particles 

Chapter X is concerned with rotational and vibrational 
contributions to the thermodyvnami« properties of gases and 
also with the subject of ortho- and pura hy lrogen The three 
following chapters present very short accounts of the statistical 
thermodynamics of radiation, adsorption and lattice defects 

Chapter XIV, on melting and crystal structure, is one 
of the best in the book and is based to a large extent on 
Prof. UsseLoupe’s own researches, The last two « hapters are 
concerned with isotopes and exchang: equilibria and with 
molecular flexibility and the statistical theory of rubber. 

There are certainly very few of the applications of statis- 
tical thermodynamics which are not in¢ luded in this interest 
ing and useful survey K. G. Dension 
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